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ABSTRACT 



The need to introduce flexibility in the optima! design of a chemical 
plant arises because of the fact that very often chemical plants have to 
operate under various different conditions. If these conditions are given as a 
discrete sequence of operating regimes, then the optimal design problem 
corresponds to a deterministic multi period problem. Also, there are situations 
when the values of some process parameters are known only approximately, 
or they actually vary or fluctuate during the operation of the plant. In such 
cases, the design strategy should be able to account for the uncertainties in 
the values of these process parameters. In both types of problems the main 
objective is to design a chemical plant that has the required flexibility to meet 
the specifications for the various conditions, while being optimal with respect 
to a selected economic index. This thesis addresses both these types of 
problems, giving the appropriate mathematical formulations and developing 
efficient solution procedures that can rigorously guarantee feasible steady 
state operation of the plant under the various conditions that may be 
encountered. 

v • 

The deterministic multi period design problem is formulated as a nonlinear 
program which has a block-diagonal structure in the constraints. One of the 
main computational difficulties faced in the numerical solution of such 
problems is the large number of decision variables involved in the 
optimization. It is shown that the existing decomposition techniques cannot be 
applied effectively to circumvent this problem. A new decomposition scheme 
based on a projection-restriction strategy is developed, which exploits the 
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mathematical structure of the problem and the fact that many inequality 
constraints become active at the optimum solution. Successful implementation 
of this strategy requires an efficient method to find an initial feasible point, 
and the extension of current equation ordering algorithms for adding 
systematically inequality constraints that become active. Systematic 
procedures are proposed to handle both these aspects. The performance of 
this projection-restriction strategy is analyzed to show that significant gains in 
computational effort can be achieved with this technique. The results of the 
numerical example show a linear increase of computational time with the 
number of periods. Thus, the above strategy provides a real possibility for 
solving the large-scale nonlinear programs that arise in the design of flexible 
chemical plants. 

In order to account for the uncertainties in the parameter values, a 
suitable strategy is to design a plant flexible enough, so that it can be 
guaranteed to satisfy the design specifications for a given range of parameter 
values, by suitable control in the operation of the plant. It is shown that a 
rigorous mathematical formulation of this design strategy corresponds to a 
two-stage nonlinear infinite program, wherein an optimization is to be 
performed on the set of design and control variables, such that the constraints 
are satisfied for every parameter value that belongs to a given bounded 
polyhedral region, while minimizing the expected value of the annual cost 
function. To guarantee feasibility of the design for every possible realization 
of the parameters, a logical constraint is incorporated in the problem. This 
feasibility constraint ensures that for every allowable value of the parameters 
(within the given range) that may be encountered during the operation, 
appropriate values of control variables can be chosen so as to satisfy the 
design specifications. In order to circumvent the problem of infinite 
dimensionality in the constraints, this feasibility constraint is shown to be 
equivalent to a max-min-max constraint, which also provides a deeper 
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understanding of the feasibility aspect of the design problem. It is shown 
that if the inequality constraint functions are convex, only the vertices of the 
polyhedral region in the parameter space need to be considered for design. 
Based on this feature, an algorithm is proposed which uses only a small 
subset of the vertices, and transforms the problem into an iterative 
multiperiod design problem. The proposed algorithm leads to an efficient 
solution procedure as shown by two example problems. 
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CHAPTER 1 
INTRODUCTION 



Chemical Plants are commonly designed for fixed nominal specifications 
such as capacity of the plant, type and quality of raw materials and products. 
Also, they are designed with a fixed set of predicted values for the 
parameters that specify the performance of the system, such as transfer 
coefficients or efficiencies and physical properties of the materials in the 
processing. However, chemical plants often operate under conditions quite 
different from those considered in the design. If a plant has to operate and 
meet the specifications at various levels of capacity, process different feeds, 
or produce several products, or alternatively when there is significant 
uncertainty in the parameter values, it is essential to take all these facts into 
account in the design. That is, a plant has to be designed flexible enough so 
as to meet the specifications even when subjected to various operating 
conditions. 

In practice, empirical overdesign factors are widely used to size 
equipment, with the hope that these factors will compensate for all the effects 
of uncertainty in the design (Rudd and Watson, 1968). However, this is clearly 
not a very rational approach to the problem, since there is no quantitative 
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justification for the use of such factors. For instance, with empirical 
overdesign it is not clear what range of specifications the overdesigned plant 
can tolerate. Also, it is not likely that the economic performance of the 
overdesigned plant will be optimum, especially if the design of the plant has 
been optimized only for the nominal conditions. 

In the context of the theory of chemical process design, the need for a 
rational method of designing flexible chemical plants stems from the fact that 
there are still substantial gaps between the designs that are obtained with 
currently available computer-aids and the designs that are actually implemented 
in practice. One of these gaps is precisely the question of introducing 
flexibility in the design of a plant. It must be realized that this is a very 
important stage in the design procedure, since its main concern is to ensure 
that the plant will be able to meet economically the specifications for a given 
range of operating conditions. Clearly, it would be desirable to consider the 
dynamic performance characteristics of the chemical processing system in such 
a design procedure so as to ensure smooth operation of the plant. However, 
in order to provide a valid framework for such considerations, it is first of all 
necessary to guarantee feasible steady state operation of the plant under the 
various conditions that may be encountered. The term flexibility here indicates 
precisely this characteristic, namely that the existence of feasible steady state 
operation is ensured for every allowable situation that is specified, and it is 
this aspect that is studied in this thesis. 

There are two classes of problems that can be considered in the design 
of flexible chemical plants. The first one is the deterministic multi period 
problem wherein the plant is designed to operate under various specified 
conditions in a sequence of time periods. Typical examples are refineries that 
handle various types of crudes, or pharmaceutical plants that produce several 
products. The second type of problem deals with the design of chemical 
plants where significant uncertainty is involved in the values of some of the 
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process parameters. Examples of this case arise when values of feed 
specifications, transfer coefficients, physical properties or cost data are not 
well established. It must be noted that in general a design problem can also 
be a combination of these two distinct types of problems. 

1.1 REVIEW OF PREVIOUS RESEARCH WORK 

Very little has been discussed in the literature about deterministic 
multiperiod problems. Loonkar and Robinson (1970), Sparrow et al. (1975), Oi 
et al. (1979), Suhami (1980) and Knopf et al. (1980) discuss design procedures 
that are applicable only to batch/semicontinuous processes. Grossmann and 
Sargent (1979) give a general formulation for designing multipurpose chemical 
plants which can also be applied to problems described by the deterministic 
multiperiod model. This involves the solution of a large nonlinear program 
(NLP), wherein the main computational difficulty is due to the large number of 
decision variables involved. Therefore, there is a clear need to develop 
efficient solution procedures for this problem. Chapters 2 and 3 present a 
decomposition scheme based on a projection-restriction strategy which is 
shown to be very efficient in solving this type of design problems. 

There have been several approaches reported in the literature on the 
problem of design under uncertainty. They differ from each other in terms of 
problem formulation as well as solution strategies, since in principle the 
problem of design under uncertainty is not well-defined. Some authors 
consider the probability distribution of the parameters as either known or 
predictable and minimize the expected value of cost. Another approach 
consists in transforming the problem to a deterministic one and assuming that 
the parameters vary within bounded ranges of values that are specified by the 
designer. A brief overview of the main features of the different approaches 
is given below. 
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Kittrel and Watson (1966) assume that probability distribution functions of 
the parameters are available, and propose to select the decision variables in 
the design so as to minimize the expected value of cost. Wen and Chang 
(1968) define the 'relative sensitivity' of the cost as the fractional change in 
the cost function from its nominal value. In selecting the optimal design they 
minimize either the expected value or the maximal probable value of this 
relative sensitivity. Weisman and Holzman (1972) incorporate a penalty in the 
cost function involving the probability of violation for individual constraints, 
and perform an unconstrained minimization of the expected value of the cost. 
Although they made an attempt to minimize the probable violation of the 
constraints, their formulation does not ensure even a lower limit on the 
probability of failure of any given constraint, which can in fact be achieved by 
using the formulation suggested by Charnes and Cooper (1959) for chance 
constrained optimization. Lashmet and Szczepanski (1974) apply Monte Carlo 
simulation for determining overdestgn factors for distillation columns. They 
perform a series of statistical experiments by choosing random values of the 
parameters (within the specified range), and in each case determine the number 
of stages in the column needed to meet the specifications. From these data 
the overdesign factor is determined as the additional number of stages 
corresponding to 90% cumulative distribution over that for nominal design. 
Freeman and Gaddy (1975) define 'dependability' as the fraction of 'time' that 
the process can meet the specifications, and use it as criterion for selecting 
the optimum design. They perform a stochastic simulation and determine the 
optimum values of the decision variables for different values of the 
parameters chosen randomly. The expected value of the cost corresponding to 
any given dependability is obtained from the results of such simulation, and 
the dependability level that minimizes the expected cost is chosen for the 
optimum design. The approach taken by Johns et al. (1976) is to select the 
optimum design by considering the effect of parameter uncertainties or time- 
varying demands on the expected value of cost computed over the whole life 
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of the plant. This procedure enables one to plan future expansions of existing 
plants and to analyze the economics of a new plant. 

It is to be noted here that in the above approaches no distinction is 
made between the two types of decision variables in the problem of optimal 
process design. The design variables, representing for example the sizes of 
equipment, get their values assigned in the design stage, and remain unaltered 
during the operation of the plant. The control variables represent the variables 
of the plant that can be adjusted in the operation. For any given design, the 
optimal plant operation itself is to be considered as a means for meeting the 
specifications while minimizing the operating cost. It requires an appropriate 
choice for the values of control variables depending on the values of the 
parameters being realized. In a realistic strategy for optimal process design 
under uncertainty, the mathematical formulation of the design problem must 
incorporate this basic difference between design and control variables. 

The following researchers have made such a distinction between design 
and control variables. Watanabe et al. (1973) apply the concept of statistical 
decision theory, by considering the problem of optimal process design as a 
two-person statistical game between Nature and the designer. They minimize a 
utility function which is a convex combination of the expected value and the 
maximum probable value of cost. That is, they follow a strategy which is 
intermediate between minimax strategy (maximal probable value of cost is 
minimized) and Bayes' strategy (expected value of cost is minimized). Avriel 
and Wilde (1969) discuss the different design strategies such as two-stage 
(here-and-now), wait-and-see and permanently-feasible programs. In a two- 
stage stochastic program, the designer selects values for the design variables 
(first stage), then observes the actual realization of the uncertain parameters, 
and accordingly chooses the appropriate values for the control variables 
(second stage). While selecting the values for design variables in the first 
stage, it is essential to ensure feasibility of the second stage sub-program, 
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namely that values of control variables can be chosen to satisfy the 
constraints. The objective is to minimize the expected value of cost while 
selecting a feasible and optimal design. Therefore, the two-stage programming 
formulation is certainly the most suitable representation for the problem of 
chemical process design under uncertainty. In the wait-and-see strategy, the 
designer waits for an observation of the uncertain parameters and then 
chooses the optimal values for both design and control variables. Here, each 
new value of the parameters results in a corresponding optimal design; or in 
other words, all decision variables are treated as control variables. In the 
permanently-feasible program, the designer selects (a single set of) values for 
both design and control variables which will be feasible for every possible 
realization of the uncertain parameters. Unlike the wait-and-see strategy, here 
the values of neither design nor control variables change with the variations in 
the values of the uncertain parameters. That is, in permanently-feasible 
program, all decision variables are treated as design variables. Avriel and 
Wilde (1969) suggest a procedure for obtaining the optimal design, which 
consists in bounding the objective function value that would be obtained at the 
solution of the two-stage program, by solving the wait-and-see program and 
permanently feasible program. However, they restrict their approach only to 
geometric programming formulations. Malik and Hughes (1980) apply a very 
similar approach for general process design problems, although there is no 
guarantee on the feasibility of the design. Also, the stochastic programming 
method they propose, based on Monte Carlo simulation, requires enormous 
computational effort. 

Takamatsu et al. (1973) assume that the parameters vary within specified 
bounds, and minimize the deviation of the objective function from its value at 
the nominal solution, while satisfying the constraints linearized around their 
nominal values. They evaluate each constraint with a separate parameter value 
that would result in the worst violation of that constraint. This formulation 
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has a major limitation apart from the fact that the constraints are 
approximated by linearizations. Their approach is to seek a design that would 
meet the specifications with a single and common operating condition for all 
parameter values. This will (if at all it does) give a very conservative design, 
probably resulting in unnecessarily expensive investment. In other words, this 
approach does not recognize and take advantage of the fact that depending on 
the values of the parameters being realized, the operation of the plant can and 
need be manipulated in order to satisfy the specifications in the most 
economical way. Dittmar and Hartmann (1976) use a similar approach as 
Takamatsu et al. (1973), but suggest instead to use the same extreme value of 
the parameter for all the constraints. They determine the design margin for 
each of these extreme values and select the largest design margin so obtained. 
Nishida et al. (1974) propose a minimax strategy wherein the maximum value 
of the cost function, as obtained at the worst parameter value in a specified 
range, is minimized by selecting the appropriate design. They view the design 
strategy as a game, wherein the uncertainty in parameter values is considered 
to be resulting in the maximization of cost, whereas the objective of the 
designer is to minimize it. It is important to note here that the design they 
come up with, actually corresponds to the optimum solution for a particular 
(namely, the economical ly worst) value of the parameters, and therefore the 
design cannot be claimed to be optimal in an overall sense when the 
parameters do take on different values. Also, the feasibility of this design 
for other values of the parameters cannot be guaranteed, since this aspect has 
not been explicitly considered in the problem formulation. Friedman and 
Reklaitis (1975a, b) deal with linear programming formulations having 
uncertainties in the coefficients. They show that this formulation can be 
applied to problems like planning future operation policies for large interacting 
systems, production scheduling, resource allocation and determining optimum 
blending schemes. They incorporate required flexibility in their system by 
allowing for possible future additive 'corrections' on the current decisions, and 
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optimize the system by applying an appropriate cost-f or-correction in the 
objective function. It is interesting to see that they were able to identify the 
need for different corrections for different outcomes of the uncertain 
coefficients, in order to make the problem feasible, and devise a procedure to 
achieve this in their computations. One obvious drawback with their approach 
is that it is applicable only to linear systems. A second limitation with their 
approach is that it is not suited for the problem of optimal design of flexible 
chemical plants, because in this case no additive corrections can be applied on 
the design variables. Instead, it is only the control variables that can be 
manipulated, so as to meet the specifications in spite of the variations in the 
values of the uncertain parameters. Kilikas and Hutchison (1980) use a linear 
process model wherein the coefficients are considered to be varying within 
specified bounds. However, their approach is not applicable for design 
purposes, because they tend to select the optimum value of their decision 
variables so as to satisfy their constraints for at least one set of parameter 
values. This objective is quite different from what is expected of a design 
strategy. In fact, they solve a steady state linearized simulation problem, with 
no mention about design variables, and therefore the optimization that they 
refer to seems to be only a suitable criterion for solution of the particular 
system of equations they present. 

With any of the above approaches, one is always faced with the question 
as to whether the designed plant can in fact be guaranteed to operate and 
satisfy the specifications for the entire range of parameter values involved. 
This question, along with the fact that the problem of optimal process design 
under uncertainty is not well-defined, requires a systematic procedure in 
formulating as well as solving such design problems. First of all, it is to be 
noted here that while minimizing the cost, the main concern of a design 
engineer must be to ensure feasible steady state operation of the plant for every 
value of the parameters within specified bounds. Grossmann and Sargent (1978) 
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propose a formulation that tries to incorporate this objective. They 
approximate the expected value of the cost by a weighted average of a finite 
number of terms, assuming discrete probabilities for a finite set of parameter 
values. They select the optimum design by minimizing this expected cost 
subject to maximizing each of the individual inequality constraints with respect 
to the parameters. In their solution procedure a small set of extreme values 
of parameters is selected by analyzing the signs of the gradients of each of 
the individual inequality constraints, and the optimization is performed for this 
set of parameter values, in the form of a multiperiod design problem. 
However, as it is shown in Chapter 5, their approach cannot always guarantee 
the feasibility requirements mentioned above. 

Among the literature on optimal design under uncertainty in other 
engineering fields, Kwak and Haug (1976) formulate the problem of 'parametric 
optimal design' wherein the current decision is chosen so as to maintain future 
feasibility even in the presence of uncertainties (with no possible corrections 
on the decision that has already been taken). This leads to a nonlinear 
program wherein the constraints correspond to individual maximization (in the 
space of the uncertain parameters) of each of the constraint functions. They 
show the application of this approach in designing structural (civil engineering) 
elements. However, the major difference here is that once the design is 
selected, no further adjustments are possible - unlike the case of chemical 
process systems. In .fact, it is interesting to see that although the freedom in 
future decisions enables a given design to remain feasible by incorporating the 
appropriate adjustments, it requires more complex mathematical representation, 
resulting in computational difficulties. In structural design it is usually the 
external load that may be uncertain, and there seems to be no adjustable 
variables available, the values of which can be chosen along with changes in 
load. The problem of design centering, tolerancing and tuning of 
electrical/electronic circuits has been studied extensively, in electrical 
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engineering literature. One objective there is to maximize the yield by 
selecting the appropriate design center and assigning the right tolerance, so 
that the effect of any uncertainties in the 'design parameters' caused during 
manufacture can be compensated by applying the required tuning. See for 
example, Bandler (1974), Bandler et al. (1975), Director and Hachtel (1977), Polak 
and Sangiovanni-Vincentelli (1979). However, there the uncertainty occurs in 
the manufacturing, and is associated with the 'design parameters' that define 
the circuits, with a difference that these can be tuned by a later finishing 
process which will try to compensate directly, for the effect of the 
uncertainties. It is to be noted here that the problem of optimal design of 
chemical plants with uncertain 'process parameters' is different in its nature 
from the above two types of design problems. This is because the 
uncertainties for this case correspond to the process parameters once a given 
design has been chosen, unlike the design of structural elements or integrated 
circuits. Hence, the methods proposed for the above types of design 
problems cannot be applied readily to the design of chemical plants. 

1.2 OBJECTIVES OF THIS RESEARCH WORK 

This thesis addresses both the deterministic mu/ti period problem and the 
problem of design under uncertainty. In Chapter 2 an efficient decomposition 
scheme based on a projection-restriction strategy is presented for solving the 
large nonlinear programming (NLP) problem that arises in formulating the 
deterministic multiperiod model. This decomposition strategy effectively 
exploits the block-diagonal structure in the constraints of this NLP, and also 
the fact that many inequality constraints become active at the solution of a 
design problem. Efficient methods have also been proposed for finding an 
initial feasible point for the design problem, and for detection of redundancy 
and selection of decision variables for the system of equations in the 
resulting optimization problem. Computational results obtained by the 
application of the proposed decomposition strategy are discussed in Chapter 3. 
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A rigorous mathematical formulation for the problem of design under 
uncertainty and its properties are presented in Chapter 4. It is shown that this 
design problem is represented by a two-stage nonlinear infinite program, 
wherein optimization is to be performed on the set of design and control 
variables so as to satisfy the constraints for ail possible values of the 
parameters (within a specified bounded polyhedral region), while minimizing the 
expected value of the annual cost function. The feasibility of the design is 
expressed by a logical constraint, which ensures that for every allowable value 
of the uncertain parameters that may be encountered during operation, 
appropriate values for control variables can be chosen to manipulate the 
operation so as to satisfy the constraints, while minimizing the operating cost. 
The main computational difficulty with this formulation arises because of the 
infinite number of constraints and infinite number of variables involved. By 
selecting a finite number of parameter values within the given bounds, and 
approximating the expected value of cost function by a weighted average, the 
number of variables can be reduced to be finite. To overcome the problem of 
an infinite number of constraints, an equivalence is established for the 
feasibility constraint which leads to a max-min-max constraint. This max-min- 
max constraint provides information regarding the extent of (in)f easibility of a 
given design. For the case of a convex feasible region, it is shown that the 
solution to the max-min-max problem lies at a vertex (corner point) of the 
polyhedron in the parameter space. Based on these features, a solution 
algorithm is presented in Chapter 5, which incorporates the projection- 
restriction strategy (as in the multiperiod design problem) in an efficient 
iterative scheme. The application of the proposed solution algorithm and the 
computational results are also discussed in Chapter 5. 
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CHAPTER 2 



THE DETERMINISTIC MULTI-PERIOD 
DESIGN PROBLEM 



There are many situations when a chemical plant has to operate under 
more than a single set of operating conditions. For example a plant may have 
to operate at various levels of capacities, process different types of feeds, or 
produce several products. This can be achieved only if the plant is designed 
with the required flexibility, so as to meet the specifications even when 
operating under various different conditions. When these conditions are either 
known or predictable, the design problem can be formulated as a deterministic 
multiperiod problem. Since in such a formulation the specifications for each 
of the various operating conditions are explicitly incorporated, the resulting 
design will naturally satisfy the feasibility requirements. The design problem 
is in fact formulated as a nonlinear programming (NLP) problem wherein the 
annual cost is minimized, subject to the constraints that represent the design 
specifications. This has been the general approach proposed by Grossmann 
and Sargent (1979), and is presented below. 
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2.1 MATHEMATICAL FORMULATION 

In the deterministic multiperiod model it is assumed that the plant is 
subjected to piecewise constant operating conditions in N successive time 
periods. Dynamics are neglected, as it is considered that the length of the 
transients is much smaller than the time periods for the successive steady 
states. The variables in this problem are partitioned into three categories. 
The vector d of design variables is associated with the sizing of the units. 
These^ variables remain fixed once the design is implemented, and do not vary 
with the changes in the operation of the plant. The vector z 1 denotes the 
control variables, the values of which can be suitably chosen in each period i 
so as to meet the specifications and also minimize the operating cost. It 
should be noted that the vector z 1 corresponds to the existing degrees of 
freedom in the operation of the plant. Finally, the vector x 1 corresponds to 
the state variables in the operating period i (i= 1,2,....N). Thus, the design 
problem leads to the nonlinear program, 



N 

minimize C = C°(d> - TcWW 

1 O M 



, z 1 ,z 2 ,...z N 
s.t. h'fdtz'tX 1 ) = 0 

i = 1,2 N 

gWrX 1 ) < 0 
1 f(d) < 0 



where d, z 1 , i=1,2,....N, are the decision variables in this problem, as the state 
variables x, i=1,2,....N, can be determined from the equality constraints which 
represent the steady state performance of the process system. Note that in 
this formulation the order in which the periods are considered can be arbitrary, 
since the operation in each period is assumed to be independent of its relative 
position in the sequence. However, any specifications involving all the periods 
can in general be represented by the last constraint, in which case f would be 
a function of d, z 1 , x 1 , i=1,2,....N. 
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2.2 COMPUTATIONAL ASPECTS 

For large industrial problems the computational requirements for solving 
the NLP in (2.1) can become rather expensive. The reason for this is that the 
number of control variables z' increases with the number of periods N, so that 
the number of decision variables in the NLP may become too large for the 
problem to be solved efficiently by the current algorithms. Since the NLP 
approach for designing flexible plants has proved to be very effective in small 
problems (Grossmann and Sargent, 1979), and it also provides a rational basis 
for overdesign, there is a very high incentive for deriving an efficient method 
for solving problem (2.1). This requires that its mathematical structure be fully 
exploited. 

In order to take advantage of the sparsity of the constraints, the state 
variables x 1 can be eliminated from the system of equations so as to reduce 
the size of the problem. This can be achieved, for instance, if the system of 
equations is ordered so as to provide a sequence of calculation where the 
number of torn variables is minimized (see e.g. Christensen, 1970). In this 
scheme, at each iteration of the optimization the ordered system of equations 
is solved. It must be pointed out, however, that by eliminating the equations 
and state variables the nonlinear program in (2.1) still has to handle the large 
number of decision variables given by d, z, i=1,2,....N. Therefore, it is 
necessary to exploit additionally another property of (2.1) for deriving an 
efficient method of solution. 

Problem (2.1) has the interesting feature that it is an NLP with block- 
diagonal structure in the constraints, as shown in Figure 2.1. Since the cost 
function is separable in the N periods this implies that if the vector d is 
fixed, then the problem decomposes in N uncoupled subproblems, each having 
as decision variables the vectors z', for i=1,2,....N. This would suggest that it 
should be possible to derive a suitable decomposition scheme which need not 
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Figure 2.1: Block-diagonal Structure in the Constraints of Problem (2.1) 
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handle simultaneously all the decision variables. Ideally, such a decomposition 
scheme should lessen the storage requirements, and more importantly it should 
reduce substantially the computational time and also enhance the reliability of 
the method for obtaining the solution. 

2.3 DECOMPOSITION STRATEGIES 

The two basic decomposition strategies that could be used for solving 
problem (2.1) are the feasible and the infeasible decomposition schemes. 

The feasible decomposition technique (Rosen and Ornea, 1963; Umeda, 
Shindo and Tazaki, 1972) consists of the following steps: 

Step 1 - Find a feasible point d, z 1 , i=1,2,....N, for problem (2.1). 
Step 2 - By keeping the vector d fixed, solve the N subproblems 



(that is, for i = 1,2,....N): 



minimize C l (d,z\x l ) 



z 



s.t. 



(2.2) 



Step 3 - Keeping the vectors z, i=1,2,....N fixed, solve the problem: 



minimize C 
d 




i=1 



s.t. 



(2.3) 



1-1,2 



g'ld.z'.x') * 0 



f(d) < 0 
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Step 4 - If convergence is not achieved, return to step 2. 

The advantage with this technique is that the original problem is replaced 
by a sequence of subproblems with a smaller number of decision variables. 
However, convergence to the solution can become extremely slow, particularly 
in the neighborhood of the solution (see Grigoriadis, 1971), since in fact this 
decomposition technique is equivalent to an alternating search in orthogonal 
directions in the space (d, (z\z 2 ,...z N ) ) 



In the infeasible decomposition technique (Brosilow and Lasdon, 1965; 
Lasdon, 1970; Stephanopoulos and Westerberg, 1975b) it is first necessary to 
reformulate problem (2.1) as: 



N 




z 1 ,z 2 z N 



s.t. 



hV.z'.x') = 



0 



(2.4) 



i = 1,2 N 



f'(d') 



< 



f(d) 



< 



c 



Since the lagrangian of this problem is given by: 



L = C°(d) ♦ ' r C'fd'.z'.x') 



i=1 
N 



+ X [ (xi)T hi + ( ^ i,T gi + (vi) 



T 



f' - ( ff ') T (d - d')] 



i=1 




(2.5) 
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where X 1 , p , v\ w , p are the Kuhn-Tucker multipliers, problem (2.4) can be 
decomposed into the following N+1 subproblems: 

minimize cW.zU 1 ) - (»') T d 1 

s.t. hV.zU 1 ) = 0 

gV.zU 1 ) < 0 
f'(d') < 0 

for i=1,2 N, and (2.6) 

N 



minimize C°(d) + ^ in') 1 d 
d i = 1 

s.t. f(d) < 0 



The infeasible decomposition strategy then consists of the following steps: 
Step 1 - Guess the multipliers r , i=1,2,....N. 
Step 2 - Solve the N+1 subproblems (2.6) 

Step 3 - If the constraints d = d', i=1,2 N, are not satisfied, adjust the 

multipliers jt' by solving the dual problem: 

N 

minimize C°(d) + V TcVd'.z'.x') + (^') T (d - d')l (2.7) 
1 2 N jft 

J7 ,77 ,..,77 1 1 



and then return to step 2. 



Note that in this decomposition technique it is not necessary to start 
with a feasible point for problem (2.1) as with the previous strategy. 
However, there are basically two difficulties when using this technique. The 
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first one is that the method may not converge to the solution due to the 
nonconvexities that are present in design problems which gives rise to dual 
gaps. This difficulty can be overcome with the method proposed by 
Stephanopoulos and Westerberg (1975a), but with the disadvantage that it 
requires a significant amount of computational effort. A further disadvantage 
with the infeasible decomposition scheme is that a feasible solution is 
obtained only at the exact solution of the dual problem. Considering that one 
is dealing with nonlinear problems, this can become a significant drawback in 
practice. 

With the two decomposition schemes that have been presented above, it 
is unclear whether problem (2.1) can be solved more efficiently than when the 
problem is tackled with all decision variables simultaneously. It is for this 
reason that an alternative decomposition strategy must be considered. 

2.4 THE PROJECTION-RESTRICTION STRATEGY 

Grigoriadis (1971) and Ritter (1973) have suggested a decomposition 
technique for mathematical programming problems with block-diagonal structure 
in the constraints, for the case when the objective function is convex and the 
constraints are linear. As per the classification given by Geoffrion (1970) this 
strategy belongs to the class of Projection-Restriction Strategies. The basic 
ideas of this strategy are described in the following steps: 

Step 1 - Find a feasible point d, z', x 1 , i=1,2,....N, for the problem (2.1) 

Step 2 - (Projection) 

Fixing the values of the design variables d, 
solve the N subproblems in (2.2) 

Step 3 - (Restriction) 

(a) For each subproblem i convert the n^ inequality constraints 
gj^ that are active in step 2 into equalities, and define 
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9 R - 9! 



i=1,2 N 



(2.8) 



where hp, g R are the redefined sets of equality and inequality 

constraints and g-j are the sets of inequality constraints 
that are not active in step 2. 



(b) Eliminate n A control variables from the vector z 



so as to define 



z R = Zj , x R 



La 



i = 1,2 N 



(2.9) 



where z R is the redefined vector of control variables which 
result from eliminating the vector z^ of n^ elements, 
and x R is the expanded vector of state variables. 



Step 4 - Solve the restricted problem: 



d,Z R ,Z R ,...Z R 



N 

minimize C = C°(d) + C'(d,z R ,x R 

i=1 



s.t. 



h R (d,z R ,x R ) = 0 



g R (d,z R ,x R ) < 0 



i=1,2 N 



(2.10) 



f(d) 



< 0 



Step 5 - Return to step 2 and iterate until 

(a) no further changes occur in the values of the variables d, or 
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(b) the same set of inequality constraints become active again, 
in step 2. 

Note that in step 4 the projection-restriction strategy really consists in 
solving problem (2.1) simultaneously for all variables, but in general with a 
smaller number of decision variables, since many of these get eliminated by 
the active constraints determined in step 2. Clearly, the effectiveness of this 
strategy relies heavily on the number of inequality constraints that actually 
become active at the solution. 

Grigoriadis (1971) and Ritter (1973) found that in their problems relatively 
few inequality constraints in the subproblems would become active. Therefore, 
they proposed to eliminate all the variables z', i=1,2,....N, in step 3, 
Grigoriadis (1971) with the use of the pseudo-inverse of the corresponding 
matrix of z, and Ritter (1973) with a square matrix which was generated when 
solving the subproblems. Unfortunately these techniques cannot be extended 
readily to the case when constraints are nonlinear, since they rely heavily on 
the assumption of linearity of the constraints. However, the basic idea of the 
projection-restriction strategy can indeed be extended to the problem of 
designing flexible chemical plants. 

An examination of the results for flexible plants obtained by Grossmann 
and Sargent (1978, 1979) shows that a surprisingly large number of inequality 
constraints are actually active at the solution, as can be seen in Table 2.1. 
The main reason for this appears to be the monotonicity of the cost function 
(in the constraint space) which seems to be characteristic of design problems. 
Since in general one can expect to have a large number of active constraints 
at the solution, it clearly suggests that the projection-restriction strategy can 
greatly simplify solving problems of the type as formulated in (2.1). However, 
for successful application of the projection-restriction strategy there are two 
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Table 2.1: Active constraints at the solution, in the design problems 
solved by Grossmann and Sargent (1978, 1979) 



Problem 



Number of 
decision 
variables 



Pipeline 



8 



Multiproduct Batch Plant 

(a) problem 1a. 10 

(b) problem 2. 14 

Reactor-Separator 

System 7 



Heat Exchange 
Network 



15 



Number of 
inequality 
constraints 



number of 
active constrai nts 
at the solution 



20 



23 
39 



11 
14 



24 



4 



65 



15 



* One of the constraints was redundant at the solution 
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problems that have to be considered. The first one is finding an initial 
feasible point in step 1. The second is a procedure for the elimination of 
variables in step 3, which avoids singularities in the system of equations. 
These are discussed in the following sections. 



2.5 FINDING A FEASIBLE POINT 

The problem of finding a feasible point for a design problem is in 
general a nontrivial task, because of the nonlinearities involved. In 
problem (2.1) the main difficulty when using the projection-restriction strategy 
consists in finding a value of d such that feasible solutions exist for the 
subproblems in (2.2). One approach to find a feasible point is to replace the 
cost function in (2.1) by the sum of squares of deviations of the violated 
constraints, thus leading to the nonlinear program: 

N M 2 
minimize 4> = ^> ^ £ max |fJ, g!(d,z',x')| 
d,z 1 ,z 2 ,...z N i=1 j«1 J 

s.t. hW.x') = 0 | (2.11) 



g , (d,z , ,x l ) < 0 
f(d) < 0 



= 1,2 N 



This problem can be handled by an NLP algorithm based on an active 
set strategy for the constraints as the one developed by Sargent and Murtagh 
(1973). Since the objective function in (2.11) has discontinuous second order 
derivatives the optimization should be performed with the steepest descent 
direction in the constraint space. As this procedure does not require an 
estimation of the inverse of the Hessian matrix, storage requirements can be 
reduced. 
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Although solving (2.11) simultaneously for all variables works very well 
for relatively small problems, it may be desirable to use a decomposition 
scheme for large problems. One alternative is to use the steps similar to 
those in the feasible decomposition strategy, with the objective function as 
given in (2.11) above. Hence it consists of the following steps: 

Step 1 - Guess a starting point d, z\x\ i = 1,2,....N. 

Step 2 - By keeping the vector d fixed, solve the N subproblems 
(that is, for i = 1,2 N): 

M 2 

minimize 

z' J=1 

s.t. hWV) = 0 (2.12) 

gW,x') < o 



Step 3 - By keeping the vectors z, i = 1,2,....N fixed, solve the problem: 

N M 2 
minimize * = 'max {o, 2 l MzW)} ] 

d i=l j=1 



s.t. 



hW.x 1 ) 



gW.x 1 ) < 0 
f(d) S 0 



(2.13) 



i=1,2 N 



Step 4 - If convergence is not achieved ($ > 0) return to step 2. 



It is observed that unlike in the case of finding the optimal solution 
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of (2.1), the convergence of this method to find a feasible point is quite good. 
The reason for this is that the NLP defined by (2.11) could have an infinite 
number of minima when the feasible region is non-empty. In this case the 
objective function defines a plateau of zero-value for the feasible region as 
shown in Figure 2.2, and hence there is usually no problem of slow 
convergence in the neighborhood of a feasible solution. Also, note in 
Figure 2.2 that outside the feasible region the contours of 4> in (2.11) are 
quadratic in the constraint functions so that the objective function will tend to 
be w^ll behaved. 



2.6 VARIABLE ELIMINATION IN THE RESTRICTION STEP 

The elimination of variables in step 3 of the projection-restriction 
strategy is performed for each period i by including the active inequality 
constraints in the set of equations. This implies that from the state and 
control variables x 1 , z 1 a new set of control variables Zp must be 

determined, and that the sequence of calculation for the new set of 
equations hp has to be derived, for each period i = 1,2,....N. 

Since a number of algorithms are available for selecting decision 
variables and determining sequences of calculation for rectangular systems 
(see e.g. Lee et al., 1966; Christensen and Rudd, 1969; Edie and Westerberg, 
1971; Leigh, 1973; Stadtherr et al., 1974; Book and Ramirez, 1976; Hernandez 
and Sargent, 1979), it would seem that they could be applied without difficulty 
in the problem of designing flexible chemical plants. It must be pointed out, 
however, that difficulties may arise when deriving the solution procedure for 
the restricted problem (2.10), since the added inequality constraints can lead to 
redundant or inconsistent equations, and hence, produce a singular system of 
equations. Therefore, these algorithms must be extended according to the 
following procedure for eliminating the variables in each i in the restricted 
problem: 
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Step 1 - Add all the active constraints g^(d,z',x') = 0, to the system of 

equations h (d,z ,x ) = 0, thus giving rise to a new system of 
equations h^d^z 1 ^ 1 ) = 0. 

Step 2 - Perform the optimal reordering of the new system of 
equations hp = 0, 

by minimizing the number of torn variables in z 1 , x 1 . 

Step 3 - Select control variables Zp as decision variables, and delete 

equations if necessary so as to obtain a non-singular square system 
of equations. 

It should be noted that due to the reordering of variables in step 2, the 

vector Zp can in fact contain some of the state variables from x. Also, it 
is essential to keep the design variables d as decision variables throughout, 
and not to force them to become either state or torn variables during the 
reordering. Again, in step 2 above, since it is possible that the resulting 
system has more equations than variables, a suitable equation ordering 
algorithm must be used, as for instance the one by Leigh (1973). The optimal 
sequence determined with such an algorithm is one where the 
system hMcLz 1 ^ 1 ) = 0, is reordered as shown in Figure 2.3 in two sets of 
equations: 

§(u,v) = 0 

(2.14) 

r(u,v) = 0. 

Here the subsystems s and r form a partition of the vector of 
equations hp, whereas the vectors of variables u and v are a partition of 

the vector [(z')"'", (x')"^]" 1 ". As shown in Figure 2.3, s is the set of non- 
recycle equations with lower-triangular structure, which can be solved 
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Figure 2.3: Structure resulting from equation ordering 
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sequentially for the vector v given a value of u, and r corresponds to the 
set of recycle equations. Since v can be treated as an implicit function 
of u, the system of equations in (2.14) can be reduced to the form: 

r(u,v(u)) = 0 (2.15) 



The vector u represents the set of decision variables Zp and torn variables 

t T , that is, u T = Efep' • * 3 and f° r the optimal sequence its dimensionality is 
at a minimum. 

4 

In order to delete the appropriate equations in step 3, it is sufficient to 
choose the largest nonsingular subset of equations at the current point. Note 
that for a fixed u, the subsystem s in (2.14) can be assumed to be 
nonsingular due to its lower-triangular structure, and the singularity of the 
system (2.14) can be analyzed through the jacobian J (r,u) of (2.15) which is 
given by 

J (r,u) = - -Sr- ) ~— (2.16 

c 3u dv x dv ' du 



This jacobian matrix can be evaluated numerically at the current point by 
performing perturbations in the vector u. To determine the equations to be 
deleted the following procedure can be followed. The square submatrix J Q of 
highest rank is obtained by performing a Gaussian elimination on the jacobian 
matrix J c in (2.16). The variables in u that correspond to the columns of the 
submatrix j£ will be chosen as torn variables t. The remaining variables in 

u will correspond to the decision variables z^. Those equations 

n 

r D in r that are not included in the rows of the submatrix J , will be 

n C 

deleted and treated as inequality constraints. In this way, the jacobian matrix 
of the resulting system of equations can be ensured to be of full rank and 
hence non-singular. Also, note that the jacobian matrix J c to be analyzed is 
usually of much smaller size than the jacobian of the system in (2.14). 
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A detailed discussion on the procedure for the selection of control 
variables in the restriction step is given in the Appendix, where it is presented 
in a more general context of selection of decision and torn variables in 
process design calculations. It is useful to note that this procedure can also 
be applied in the initial stage of formulating the design problem (2.1), in order 
to identify the proper choice for the decision variables in the design problem 
and also to detect any redundancy or inconsistency in the system. This is an 
important aspect in formulation of problems for the design of flexible 
chemical plants since one has to be careful to maintain a positive number of 
degrees of freedom, and also have a consistent system. 

The procedure indicated above for finding an initial feasible point and for 
the variable elimination in the restriction step, complete the algorithm required 
for the projection-restriction strategy. The next chapter presents a discussion 
on the expected performance of the proposed decomposition strategy for 
solving deterministic multiperiod problems, and an example problem which 
illustrates the application and confirms the expected performance of the 
projection-restriction strategy. 
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CHAPTER 3 



NUMERICAL PERFORMANCE OF 
THE PROJECTION-RESTRICTION STRATEGY 



The projection-restriction strategy presented in the previous chapter is 
expected to be very efficient, provided that a large fraction of the original 
inequality constraints actually become active at the solution of the design 
problem. This is in fact true in most of the design problems, as shown in 
Table 2.1, thus making the above strategy an effective decomposition scheme 
in solving large multiperiod problems that arise in the design of flexible 
chemical plants. This chapter presents an analysis of the relative gain in 
computational effort obtained by using the proposed decomposition strategy, 
and the effect of the size of the design problem as well as the number of 
active constraints on the expected numerical performance. An example 
problem has also been solved to illustrate the application of the projection- 
restriction strategy, and to observe the actual performance. 
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3.1 ANALYSIS OF THE EXPECTED NUMERICAL PERFORMANCE 

The relative advantage in using a decomposition technique may in general 
depend on the size of the problem being solved, and hence it is appropriate to 
analyze the characteristics of the proposed decomposition strategy from this 
viewpoint. 

Based on some simple assumptions, a relationship can be derived 
between the CPU-time and the size of the design problem in terms of number 
of periods. 

Let n^ be the number of design variables, n z the number of control 
variables in each period and N the number of periods considered in the 
design problem. Assume that the CPU-time t p for the projection step is 
given by 



t p = a p 



(n z ) P (3.1) 



If a, 0 ^ a ^ 1, is the average fraction of the control variables remaining in 
the restriction problem, then the CPU-time tp for the restriction step can be 
expressed as 



t R = a R (n d ♦ N a nj' 



(3.2) 



Also, let the CPU-time tQ for solving the problem without decomposition be 

l Q = a Q ( n d + N ^ (3 - 3) 



If K is the number of iterations (passes) through the projection (and 
restriction) steps that is required for convergence, then the total CPU-time 
needed to solve the design problem using the decomposition strategy is 



t pR = K N a p (n z ) P * K a R (n d ♦ N a nj' 



(3.4) 
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In general the exact values of p, q, r and a p , aQ, a p , K depend on the 
particular problem at hand. However, for a given problem the values 
of ap, aQ, a p can be expected to be of the same order of magnitude, and 
for a gradient based non-linear programming algorithm one can expect to have 
the values of p, q and r to lie between 2 and 3. Also, the value of K is 
likely to be small. 

If all the control variables are eliminated in the restriction step, a = 0 
and from (3.4) it is clear that the CPU-time t pp is linear in N the number 
of periods, as given in (3.5) below: 

a = 0: t pR = K N a p (nJ P + K a R (n^' (3.5) 

For a > 0 the relative savings in computational time in using the 
decomposition strategy can be determined from (3.3) and (3.4), 

t pR K N a p (n 2 ) P + K a R (n d + N a 

- = - — (3.6) 



J Q ( n d + Nn z) 



Since a p , a R , aQ are of the same order of magnitude; 0 < a < 1, 
and p, q, r > 1, it is clear from (3.6) that the relative advantage in using the 
decomposition strategy is enhanced by larger values of N and smaller values 
of a. In fact, for a given value of N there is a threshold 
value a X , below which savings in CPU-time can be ensured by using the 
decomposition strategy. This threshold value determines a useful range 
0 < a < a*, which can be determined from (3.6) with t pp ^ tQ, thus 

obtaining a X as 

«' ■ ^[{(^)(^)M^)(n Z )V'-"J 



Assuming a R = a R = 8q and p = r = q, (3.7) can be simplified as 
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Number of Periods, N 



Figure 3.1: Minimum fraction of control variables that must be 
eliminated as given in (3.8) 
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(3.8) 



where j3 indicates the minimum fraction of control variables to be eliminated 
in the restriction problem. Figure 3.1 shows some plots of /?* versus N as 
given by the expression in (3.8) above for the case when n^ = n , There are 
two sets of three plots, for q = 2, 3 and K = 1, 2, 3. As can be seen from 
this Figure 3.1, fi is smaller for larger values of q and for smaller values of 
K. Also, for a given q and K, /? t decreases rapidly with N even for 
relatively small values of N, and approaches zero asymptotically for large N. 

As an example, take the case when q = 2 and when a single pass 
(K = 1) in the decomposition strategy can lead to the optimal solution. If the 
design problem involves five periods, only 8% of the control variables must 
be eliminated to achieve a relative gain in CPU-time with the decomposition 
strategy. These percentages increase to 18% for K = 2 and 28% for K = 3. 
For a ten period problem the percentages reduce respectively to 4%, 10% and 
15% for K = 1,2 and 3. Thus, in general one can expect to obtain significant 
gains in computational time when solving multiperiod design problems with the 
projection-restriction strategy, even when the number of active constraints is 
not very large. 

3.2 NUMERICAL EXAMPLE 

To confirm the expected performance of the proposed projection- 
restriction strategy, an example problem is solved. The flowsheet consists of 
a reactor and a heat exchanger as shown in Figure 3.2. The reaction is 

assumed to be first order exothermic, of the type A -> B. The flowrate 
through the heat exchanger loop is adjusted to maintain the reactor temperature 
below T' lmax as given in Table 3.1 and to get a minimum 
of 90% conversion. 
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This plant is to be designed so as to produce different products 
N successive periods within each year. The performance equations of 
such a system, for any period i, i = 1, 2, ... N, are as follows: 



in 



Reactor, material balance: 



F o < C Ao * C A1 » ' C Ao = V k o expt-E'/RTVC^, 0.9) 
Reactor, heat balance: 

<- A <n F o ( C Ao * C A1 > ' C Ao = F o C p < T 1 - V + Q HE < 3 - l0 > 



Heat exchanger, heat balances: 

Q HE = F 1 C p (T 1 " T 2> (3 - 11> 
Q HE " F w C pw (T w2 " T w1> (3 - 12) 



Heat exchanger, design equations: 

Clj_l£ = A U (AT)^ (3.13) 
(AT ,i ( T 1 ~ T w2) '( T 2" T wl) 



The data for the problem are given in the Table 3.1. The design problem 
corresponds to an optimization problem with the equality constraints given by 
the performance equations above, and the following inequality constraints for 
each period i = 1, 2, ... N. 

V > 0 (3.15) 
A > 0 (3.16) 
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V 1 > 0 (3.17) 



V - V 1 > 0 (3.18) 

"4 * 0 (3-19) 

f!| > 0 (3.20) 

°- 9 * < C Ao - C A1>/ C Ao * 1 -° (3 - 21) 

T 1 S T 1max < 3 - 22 > 

t!, - > 0 (3.23) 

T w1 5 T w2 * 356 (3 ' 24> 

T 1 , - tL, ^ 11.1 (3.25) 



1 w2 

' - T 

'2 ' w1 



To " T.„, > 11.1 (3.26) 



The objective function being minimized is the total annual cost ($/yr), 

N 

C = 691.2 V 0 - 7 ♦ 873.6 A 0,6 + jT «' (l.76 ♦ 7.056 f!,) (3.27) 

i=1 



where t' corresponds to the number of hours of operation for each period 
i=1,2,....N, in one year.' The objective function includes the investment cost of 
the reactor and heat exchanger, and the operating cost of the cooling water 
and recycle. 

There are 2 + 9N variables V, A, C^ r T' r tL tJ^,, F' v F^, V' ( (AT)L 
Q^lg, i = 1,2,....N; 6N equations and 2 + 10N inequality constraints and 
bounds, for a problem with N different periods. This gives rise 
to 2 + 3N degrees of freedom for the design problem. Selecting as decision 
variables the design variables V, A, and the control variables TL T^- 
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i=1,2,....N, the sequence of calculation for the equations in each period is 
given in Table 3.2. The corresponding nonlinear program consists of 2 + 3N 
decision variables, 6N nonlinear inequality constraints and N linear inequality 
constraints. Note that (3.15), (3.16), (3.22), (3.24) and (3.26) are simple bounds 
on the decision variables, that (3.23) is a linear inequality and the remaining 
constraints are nonlinear. 

The problem has been solved for five cases corresponding 
to N=1,2,3,4,5. In each case, the plant is designed to produce N different 
products that have different feed rates, concentrations, reaction rate constants, 
etc. as indicated in Table 3.1. In all the five cases, when solving the 
subproblems in the projection step it is found that constraint (3.21) is active at 
its lower bound, and constraints (3.22) and (3.24) are active at their upper 
bounds, for all periods. Adding these active constraints to the equations in 
Table 3.2, the variables T 1 -, Ti->, i=1,2,...N were eliminated by ordering 

the new system of equations. This gives rise to only two decision variables 
V and A for the restricted problem, as shown in Table 3.3. 

The starting point given in Table 3.1, which is infeasible, is used for all 

five cases. The initial feasible points used in the projection-restriction 

strategy were obtained by minimizing alternately with respect 

to d and z, i=1,2,....N, the sum of squares of deviations of violated 

constraints. The optimizations were performed using the variable metric 

projection method (Sargent and Murtagh, 1973), and the solutions were obtained 

-2 

with a tolerance of 10 for the norm of the gradient of the objective function 
projected in the constraint space. 

The optimal sizes of the reactor and the heat exchanger are presented in 
Table 3.4 for the five cases. The formulation of the problem itself ensures 
that these optimal designs are flexible, as they meet the specifications for the 
various products involved at a minimum annual cost. 



- KpH - 



- 39 - 



The computing requirements for finding the initial feasible points are 
shown in Table 3.5. Here, it was found that optimizing alternately 
for d and z 1 is more efficient than considering all these variables 
simultaneously, particularly when the number of periods is large. However, the 
more significant gains in computational requirements are achieved when the 
projection-restriction strategy is applied, once the problem becomes feasible. 
Table 3.6 gives the CPU-time requirements for the projection and restriction 
steps. Table 3.7 and Figure 3.3 give a comparison of the computational 
requirements in solving the design problem with and without the use of the 
decomposition scheme. A striking feature in the performance of the proposed 
decomposition strategy is that the CPU-time increases only linearly with the 
size of the design problem. It is interesting to note that the design problem 
for the five-period case was solved by using the proposed decomposition 
strategy in only 31.4 sec which is about the same time required for the one- 
period problem without using any decomposition. 

3.3 DISCUSSION 

The results of the above example show that the performance of the 
proposed decomposition strategy for the design of flexible chemical plants is 
very encouraging. In fact, the decomposition strategy has been applied to 
solve the design problems in Chapter 5 as well, requiring only moderate 
amount of computational time. It is to be noted in Table 3.7, that the 
computational effort required to solve the design problem of this chapter has 
been reduced by an order of magnitude with the application of the projection- 
restriction strategy. An important trend in the results is that the computational 
time required is approximately linear with the size of the design problem 
(number of periods). This suggests that the reduction in computational effort 
with the proposed method is even more dramatic in larger problems. This is 
to be expected, since experience with different nonlinear programming 
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algorithms indicate that they are much more likely to be successful in 
converging to the optimal solution when the number of decision variables is 
relatively small. Also, it will lead to a more reliable method in obtaining the 
desired optimum design solutions, since the large nonlinear programming 
problem is reduced to a sequence of smaller subproblems. 

Note that the proposed decomposition strategy is conceptually quite 
general, and does not presuppose the use of any particular optimization 
algorithm for solving the resulting nonlinear programming (NLP) problems. 
Although no numerical experimentation has been performed with different NLP 
algorithms, it is clear from the derivations given at the beginning of this 
chapter, that one can expect significant savings in computational effort with 
almost any algorithm for solving the nonlinear program. The exponential 
increase in CPU-time with the increase in number of periods that is observed 
for solving the example problem without decomposition is typical of any 
optimization algorithm. In fact, the above example problem has been solved 
for the case of 1, 2 and 3 periods (without the application of any 
decomposition strategy) using the Han-Powell Algorithm (Powell, 1977) for 
optimization; requiring 14.1sec, 126.5sec and 543.3sec respectively to attain the 
solution (as against 32.2sec, 176.2sec and 479.6sec respectively when using the 
variable metric projection algorithm). These computational results seem to 
compare with those obtained by using the variable metric projection algorithm 
of Sargent and Murtagh (1973), although in general the latter was more 
accurate because of its stringent termination criteria - for any given error 
tolerance. From these results it can be seen that whatever be the optimization 
algorithm used, the application of the decomposition strategy will certainly 
enhance the efficiency in solving the multiperiod design problems. 
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Figure 3.2: Flowsheet for the Example Problem 
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Table 3.1: Data for the Example Problem 



Period i: 


1 


2 


3 


4 


5 




(E'/R) 


555.6 


583.3 


611.1 


527.8 


500.0 


K 


<- AH >rxn 


23260 


25581 


27907 


20930 


18604 


kJ/kgmole 


K o 


0.6242 


0.6867 


0.7491 


0.5619 


0.4994 


hr" 1 


c p 


167.4 


188.4 


209.3 


146.5 


125.6 


kJ/kgmole K 




32.04 


40.05 


48.06 


24.03 


32.04 


3 

kgmole/m 




45.36 


40.82 


36.29 


49.90 


54.43 


kamole/hr 


r\ 

1 max 


389.0 


383.0 


378.0 


394.0 


400.0 


K 


U 


1635.34 


1635.34 


1635.34 


1635.34 


1635.34 


kJ/m z hr K 


' o 


333.3 


333.3 


333.3 


333.3 


333.3 


K 




300.0 


O \J \J . u 


JUU.U 






K 

In 


Starting Point: 


V 


B 


14.16 


m 3 








A 




1 1.15 


m 2 












367 


K 












328 


K 


i-1,2 N 






' w2 




333 


K 
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Table 3.2: Equation Ordering in the Projection Step 



VARIABLES 





























V 


A 


A 




'w2 


< AT >m 


Q HE 


c' 

^A1 


V' 




F w 


14 






X 


X 


X 


X 












13 




X 








X 


X 










10 






X 








X 


X 








S 






X 










X 


X 






11 






X 


X 






X 






X 




12 










X 




X 








X 
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Table 3.3: Equation Ordering in the Restriction Step 



VARIABLES 





\i 

V 


A 


c' 

L A1 


1 1 


1 w2 


V 1 

V 


q' 

HE 


(AT)' T' 
(AU m '2 


F i 

1 


F i 

h w 


21 






X 
















22 








X 














24 










X 












9 






X 


X 




X 










10 






X 


X 






X 








13 




X 










X 


X 






14 








X 


X 






X X 






1 1 








X 






X 


X 


X 




12 










X 




X 






X 
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Table 3.4: Solution of the Example Problem 



* 

Number of 
Periods N 


V m 3 


A m 2 


Optimum Annual 
Cost $/yr 


1 


5.318 


7.562 


0.980 X 10 4 


2 


5.318 


8.417 


1.010 X 10 4 


3 


5.318 


9.513 


1.042 X 10 4 




7.915 


9.262 


1.096 X 10 4 


5 


7.915 


9.095 


1.080 X 10 4 



* Number N indicates periods 1,2,....N taken together. 
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Table 3.5: Computational results for finding an initial feasible point, 

CPU-Time (DEC-20): 



* 

Number of Optimize all Variables Optimize Alternatively 

Periods N Simultaneously d and z 



1 0.228 0.204 

2 0.371 0.432 

3 3.745 0.670 

4 1.763 0.878 

5 2.954 1.083 



* Number N indicates periods 1,2,....N taken together. 
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Table 3.6: Computational results for the Projection-Restriction Algorithm, 

CPU-Time (DEC-20): 



* 

Number of 
Periods N 


Projection 


Restriction 


Total CPU-Time 
DEC-20 sec 


1 


4.531 


1.304 


5.835 


2 


12.129 


2.221 


14.350 


3 


15.918 


1.416 


17.334 


4 


21.998 


2.181 


24.179 


5 


27.311 


3.005 


30.316 



* Number N indicates periods 1,2,....N taken together. 
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Table 3.7: Computational results for solving the design problem, 

CPU-Time (DEC-20): 



* 

Number of Number of Number of Computational Time 



Periods N Decision Inequality without with 

Variables Constraints Decomposition 

1 5 7 32.2 6.0 

2 8 14 176.2 14.8 

3 11 21 479.6 18.0 

4 14 28 — 25.0 

5 17 35 — 31.4 



# Number N indicates periods 1,2,....N taken together. 
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Figure 3.3: Computational time for solving the design problem 
versus the size of the problem 
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3.4 NOMENCLATURE FOR THE NUMERICAL EXAMPLE PROBLEM 



Heat transfer area of the heat exchanger, m 

3 

Concentration of reactant in the feed stream, kgmole/m 

Concentration of reactant in the product stream, kgmole/rrT 
Heat capacity of reaction mixture, kJ/kgmole K 



Ratio of activation energy to gas constant, K 
Feed flowrate, kgmole/hr 

Flowrate of the recycle, kgmole/hr 
Flowrate of the cooling water, kg/hr 

Heat of reaction, kJ/kgmole 

Arrhenius rate constant for reaction, hr 
Total number of periods considered for design, (integer) 
Heat exchanger load, kJ/hr 
Temperature of feed stream, K 

Reaction temperature, K 

Recycle temperature, K 
Inlet temperature of cooling water, K 
Outlet temperature of cooling water, K 

Mean temperature difference, K 

o 

Overall heat transfer coefficient, kJ/m hr K 

3 

Reactor volume (design capacity), m 

3 

Reaction volume in the operating period i, m 

weighting factor for the period i. 
Superscript i refers to the period of operation, i=1,2,....N 
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CHAPTER 4 

OPTIMAL PROCESS DESIGN UNDER UNCERTAINTY 



In the optimal design of chemical processes it is very often the case 
that considerable uncertainty exists in the value of some of the parameters. 
For instance, these parameters may correspond to internal process parameters 
such as transfer coefficients, reaction constants, efficiencies or physical 
properties. In addition, the uncertain parameters may also include those 
external to the process itself, such as specifications in the feed streams, 
utility streams, environmental conditions or cost data. The values of these 
parameters are either known only approximately during the design stage, or 
they may actually vary or fluctuate during the operation of the plant. Hence, 
in designing an optimal chemical process system, these parameters have to be 
treated as having uncertain values, and the effect of these uncertainties on 
both the optimal ity and the feasibility of operation must be considered. 

To account for the uncertainties in the values of these parameters, 
normally the procedure that is used in practice is to assume some nominal 
values for the parameters, and then apply empirical overdesign factors to the 
resulting design. Since this procedure lacks a firm rational basis, a number of 
different methods have been suggested to account for the uncertainties in a 
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more systematic manner. These methods differ mainly in the basic design 
strategies that are postulated, since in principle the problem of design under 
uncertainty is not well-defined. A brief review of the different approaches 
and their shortcomings is given in Chapter 1. 

As mentioned in Chapter 1, the main objective in the problem of optimal 
process design under uncertainty is to guarantee feasible steady state operation of 
the plant under the various different conditions that may be encountered, so as to 
meet the specifications. In other words, while minimizing the cost to select the 
optimal design, the main concern of a design engineer must be to ensure that 
even when the parameters take on different values within a given bounded 
region, suitable manipulation in the operation of the plant can be achieved in 
order to satisfy the specifications. Grossmann and Sargent (1978) have 
proposed a design strategy which tries to incorporate this objective by 
designing a plant flexible enough to accommodate the uncertainty in the values 
of the parameters. The basic idea in this strategy is to take advantage of the 
fact that control variables can be adjusted to satisfy the design specifications 
during the operation of the plant, as it is only the design of the plant itself 
that will remain fixed. Based on this strategy, a rigorous mathematical 
formulation of the problem of optimal design of flexible chemical plants with 
uncertainty in parameter values is given in this chapter. As will be shown in 
Chapter 5, this formulation is more general than the one presented by 
Grossmann and Sargent (1978) which can fail in some cases to guarantee 
feasibility. The interpretations and properties of the proposed formulation for 
the case of nonlinear convex functions are also presented here, whereas the 
proposed solution algorithms with numerical examples are given in Chapter 5. 
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4.1 MATHEMATICAL FORMULATION 

The variables in the design problem of a chemical plant with uncertainty 
in parameter values can be partitioned into four categories. The vector d of 
design variables is associated with the sizing of the units. These remain fixed 
once the design is implemented, and do not vary with the changes in the 
operation of the plant. The vector z denotes the control variables, the values 
of which can be chosen during the operation so as to meet the specifications 
and also to minimize the operating cost. The vector x corresponds to the 
state variables which are determined by solving the set of equations 
representing the process system. Finally, 8 is the vector of independent 
parameters in the design whose values are subject to uncertainty. Assuming 
that bounded range of values are specified for these parameters, the region T 
that is defined to contain all possible values of these parameters is given by 



where 9 and 9 represent given lower and upper bounds on 8. We assume 
that the control variables z lie in a compact set. 

In order to derive the mathematical formulation, it is convenient to 
consider the design strategy as being composed of two stages: an operating 
stage and a design stage. 

(1) Operating stage: Assuming that a given design d has been selected, it is 
anticipated that the plant will be operated optimally while satisfying the 
constraints of the process for all possible realizations of the parameters 8 
within the region T. Hence, the objective in this stage is to select for every 

realization 8 6 T, a control z which is optimal and feasible. 

Clearly, for the given d and for any value of 8, the state variables can 
be expressed as an implicit function of the control variables z, from the 
system of equations of the process, 



T 




(4.1) 
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h(d,z,x,0) = 0 x = x(d,z,0> (4.2) 

Since the control variable z should be selected so as to satisfy the 
specifications given by the vector of inequality constraints, 

g(d,z,x,0) = g(d,z,x(d,z,0),0) = f(d,z,0) < 0 (4.3) 

the optimal feasible operation of the plant that minimizes the cost will be 
given by the nonlinear program 

minimize C(d,z,0) 
z 

s.t. f(d,z,0) < 0 (4.4) 



Since the control variables z lie in a compact set, the minimum is attained. 
The solution of this problem defines a cost function C*(d,#) which corresponds 
to the optimal operation of the plant for fixed values of d, 6. Moreover, if 

the optimization is performed for every realization 6 £ T, the average cost of 
operation will be given by the expected value £ {c*(d,0)j. 

(2) Design stage: In order to achieve the basic objective of feasibility of 
operation in the region T of parameters, the design variable d must be chosen 

so as to ensure that for every value 0 £ T the control variable z can be 
selected to satisfy the constraints in (4.4). Note that an improper selection of 
d can lead to infeasible operation for some realizations of 6 in which case no 
selection of the control z will exist so as to satisfy the inequality constraints 
in (4.4). Furthermore, in order to achieve the optimal design, the design d 
must be selected so as to minimize the expected value of the optimal cost 
function C {6,6) over the entire region T. 
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This strategy can then be expressed mathematically as 

minimize E { C*(d,0) } 
d 8ET 

s.t. V^^T{3z(Vjej[ f-{d.z.O) < 0 ])} (4.5) 

where J = {1,2 m} is the index set for the components of vector f. The 

constraint in (4.5) is denoted as the feasibility constraint, because the existence 

of a feasibility region of operation for every #ET can be ensured if and only 
if this constraint is satisfied. In fact, this logical constraint states that for 

every point 6 £ T in the space of parameters, there must exist at least one 
value for the vector z of control variables that gives rise to non-positive 
values for all the individual constraint functions. Qualitatively, this means that 
irrespective of the actual values taken by the parameters, with the selected 
design d, the plant can be operated to satisfy the specifications. 

Since the objective function in (4.5) is itself determined through the NLP 
in (4.4), the problem of optimal design under uncertainty can be formulated in 
its final form as a two-stage programming problem: 

minimize E { min C(d,z.0) | f(d,z,(9) < 0 } 
d 6E1 z 

s.t. V^ T {3 2 (Vi GJ [ f j< d ' z '^ - 0 ])} (4 - 6) 

Note that since there are an infinite number of possible realizations for 
the values of the parameters 6, and since the optimal operation of the plant is 
implicitly dependent on 6, the overall number of decision variables involved in 
problem (4.6) is infinite. This is because for every value of 0 an optimal 
value of the control variables z is being chosen. Also, note that the 
feasibility constraint represents an infinite set of constraints since the 
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inequalities in (4.4) are defined for the infinite set of values SET. Therefore, 
problem (4.6) corresponds to a two-stage nonlinear infinite program (NLIP). 

4.2 SIMPLIFICATION OF THE TWO-STAGE NLIP 

The two-stage nonlinear infinite program in (4.6) that represents the 
mathematical formulation of the design problem under uncertainty poses great 
computational difficulties for numerical solution, and in fact has a more 
complex structure than the semi-infinite programs treated in the literature (see 
Hettich, 1979). A first step in simplification so as to make the problem more 
tractable, is to perform a discretization over the parameter space in order to 
approximate the expected cost by a weighted average (Grossmann and Sargent, 
1978), which reduces (4.6) to: 

N 

minimize y w { C(d,z 1 ,^ 1 ) 
d,z 1 ,z 2 z N PT 

s.t. f(d,z\$ l ) £ 0, i = 1,2 N (4.7) 

V^i{3z(Vi€j[ fj(d,z,<?) < 0 ])} 



where the weights a;' correspond to discrete probabilities for the selected 

finite number of parameter points £ T, i=1,2,....N. Since very often the 
probability distribution functions for the parameters are not available at the 
design stage, these weights may be selected by the designer so as to reflect 
subjective probabilities. 

With the above simplification, the number of decision variables in (4.7) is 
finite, since optimization is performed over the vector d of design variables 

and the finite number of vectors z\z^, z^ of control variables. The control 

variables z' are selected to satisfy the corresponding constraints f(d,z 1 ,^ 1 ) £ 0, 
so as to result in an optimal feasible operation at the point d x of the 
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parameter space. 



Note that despite this discretization, the feasibility 



constraint is still imposed so as to restrict the choice on the design d to 

guarantee feasible operation for every point 6 G T. Thus, the formulation in 
(4.7) is a nonlinear semi-infinite program (NLSIP), with a finite number of 
decision variables and infinite number of constraints. 

It is interesting to note that if the feasibility constraint is excluded in 
(4.7), the resulting structure of the problem is equivalent to that of a 
deterministic multiperiod problem, where the plant operates in each period with 
the parameter value and with length of period proportional to w\ As 
discussed in Chapters 2 and 3, this class of multiperiod design problems can 
be solved very efficiently with the projection-restriction strategy. The 
question that immediately arises then is whether a finite number of points in 
#-space can be selected so that by ensuring feasibility of the design for those 
points, one can guarantee that the feasibility constraint in (4.7) will be 
satisfied. If such a choice of finite number of parameter values were possible 
one could clearly solve problem (4.7) as an equivalent deterministic multiperiod 
design problem. In order to answer this question, it is essential first to 
reformulate the feasibility constraint in (4.7) so as to make it more amenable 
for analysis. 

4.3 REFORMULATION OF THE FEASIBILITY CONSTRAINT 

The logical constraint 



which ensures overall feasibility of the design is the major source of 
computational difficulty in numerical solution of the problem represented by 
the NLSIP in (4.7). The reason is that it involves an infinite number of 




(4.8) 
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inequality constraints for which feasibility has to be tested. The following 
theorem provides a possibility for circumventing this problem. 

Theorem 1. The logical constraint (4.8) and the max-min-max constraint, 



max min max f.(d,z,0) ^ 0 (4.9) 
del z jGJ J 



are exactly equivalent. 



Proof: This theorem can be proved in two parts, namely, 

(4.8) => (4.9) and (4.9) (4.8), 
as given by Polak and Sangiovanni-Vincentelli (1979), and by Halemane and 
Grossmann (1981). An alternative proof which is simpler and more direct is 
given here. By the definition of the terms and relationships used in (4.8) and 
(4.9), the following equivalences apply, by considering global max and min 
operators: 

V^ T {3 2 (Vj€j[ fj<d,z,0> < 0 ])} 
<==> Y ^€T { 3 z( max f.(d,z,0) < 0 )} 

Y 6£1 { min max fAd,Z.6) < 0 } 

z jGJ 

£z^> max min max f.(d,z,0) ^ 0 

06T z jGJ J 



From these steps the equivalence of the first and last relations is established. 

Q.E.D. 
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With this alternative and equivalent formulation of the feasibility 
constraint the optimal design problem in (4.7) can be rewritten as 

N 

minimize y w ] C(d,z\d*) 

d, z \z 2 z N pr 

s.t. f(6,z\9 l ) < 0, i = 1,2 N (4.10) 

max min max f-(d,z,0) ^ 0 
8ET z jGJ J 

In order to describe qualitatively the significance of the max-min-max 
constraint, note that the inequality constraints are satisfied for non-positive 

values of the constraint functions fj(d,z,#), j 6 J. Hence, the worst constraint 
function is that which is most likely (to be) violated, and is denoted by the 
index j which corresponds to the maximum valued function fr(d,z,0) for given 

d, z, 8. The control z that minimizes this function Md,I,0) for any given d, d, 
corresponds to the most feasible operation (even) for the worst constraint. 
Then, a critical parameter value 6° can be defined as that for which the worst 
constraint fj is maximized while having the control z, for a given design d. If 

for a design d, a control variable z can be chosen to satisfy the constraints at 
every critical parameter value 6°, then the design d can be guaranteed to have 

feasible operation at every 6 E J. 

Note that in the above formulation the max-min-max constraint provides 
the possibility of circumventing the problem of simultaneously handling the 
infinite number of inequality constraints. The reason is that the max-min-max 
constraint determines a point 6° for which the inequalities are most likely to 
be violated, while requiring that these inequalities be satisfied at that point. 
However, this constraint involves solving the max-min-max problem, 

max min max f.(d,z,0) (4.11) 
del z jGJ J 
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which is structurally far more complex than the min-max problems generally 
addressed in literature (Danskin, 1967; Demyanov and Malozemov, 1974). The 
main computational difficulty arises because of the non-differentiability of the 
max-min-max functions as well as the difficulties associated with evaluation of 
such complex functions (Poiak and Sangiovanni-Vincentelli, 1979). Furthermore, 
it is not clear under which circumstances the solution of this subproblem is 
unique, since in Theorem 1 global max and min operators had to be used for 
the proof. Therfore, it is desirable as a next step to examine the properties 
and interpretation of the max-min-max constraint so as to gain a better insight 
and understanding from it. 

4.4 INTERPRETATION AND PROPERTIES OF THE MAX-MIN-MAX CONSTRAINT 

The max-min-max constraint can be written as a constrained max-min 
problem, by introducing an extra variable u, 

max min { u | u > f.(d,z,0), VjGJ } ^ 0 (4.12) 

eeJ z J 

It then follows that for a given point SET the value f(d,0) determined by 

f{d,0) = min { u | u > f.(d,z,0), Vi^ J } (4 - 13) 
z 

indicates the extent of (in)feasibility of operation of the design d for that 
point 6. A negative value of f{6,6) indicates a finite region of feasibility and 
a positive value indicates inf easibility. Thus, the value f(d,0) can be 
interpreted to be a good measure of ( i n)feasi bi I ity of operation at the chosen 

point 6 £ T. Since the constraint in (4.12) leads to a point 6 C which maximizes 
f(6,6), 6 C corresponds to a critical point in the parameter space for which the 
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design d has either the smallest degree of feasibility (if f{6,6 ) ^ 0), or the 
largest degree of inf easibility (if f(d,6 c ) > 0). 

To illustrate these ideas consider the following set of two constraints 
which involve one variable z and one parameter 8: 

f 1 * -z + e < o 

f 2 = z - 20 - d + 2 < 0 (4.14) 
1 < 6 < 2 

*• 

Figure 4.1a shows a plot of the feasible region on 6-z space, for a design 
corresponding to d = 0.5. As can be observed in the figure, the size of the 
feasible region increases as 6 increases, with 6 = 1 being infeasible, 8 = 1.5 
being feasible at one single value of z, and 0 = 2 having a finite region of 
feasibility. This can also be observed from Figure 4.1b where the feasible 
region on z-u space is shown for these three different 0-values. The value 
of f is determined by solving the problem, 

f(d,8) = min u 
z 

s.t. f 1 « -2 ♦ 6 < u (4.15) 

f 2 = z - 26 - d + 2 < u 
1 < 6 < 2 

for d = 0.5 and 1 < ^ < 2, and its results are plotted in Figure 4.1c. Note 
that f = 0 for 6 = 1.5 which has a single point of feasibility. Also, negative 
values of f correspond to finite regions of feasibility as for instance at 
6 = 2, and positive values of f are associated with inf easibility as for 
6 = 1.0, which is a critical point where the maximum of f is attained. Note 
also that f decreases monotonically with increasing 6, while the feasible 
region gets expanded. From these observations it is clear that for a given 
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Figure 4.1: Feasible region and fid,0) for constraints (4.14) with d=0.5 



f 





JL V 



\ 



(c) 



e = 1 




e = 1.5 



o 









1 \- mm 


A 




l >v 






Z 









6 = 2 




f 7> 



(b) 



Figure 4.2: Feasible region and p<d,0) for constraints (4.14) with d=1.0 
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design d and for any chosen parameter value 6, the value of -f(d,6) can be 
interpreted as a measure of the depth of the feasible region, in the z-u space, 
for operation. This region corresponds to the projection of the actual overall 
feasible region in the d-0-z-u space onto the z-u space for fixed values of d 
and <9. 

To study the effect of changes in d f the region of feasibility for (4.14) is 
shown in Figure 4.2a and 4.2b for d ■ 1.0, and the corresponding plot of f is 
given in Figure 4.2c. By changing d from 0.5 to 1.0, overall feasibility has 
been achieved for all values of 6 in the specified range 1 < 0 < 2. Again, it 
is clear from Figure 4.2a and 4.2b that 6 = 1 is a critical point, since it 
corresponds to the maximum value of y-(d,#) for the specified range 
1 < # < 2. Thus, the design d = 1.0, which is feasible for the critical point 
6 = 1 is found to be also feasible for the entire range 1 < 0 < 2. 

The example above would suggest that feasible operation in the design 
can be guaranteed by considering one single critical #-point. The existence of 
a single critical parameter value can be associated with the monotonicity of 
f(d,9) with respect to Q, which is in fact observed in this example (see the 
plot of f versus 8 in Figure 4.2c). However, this may not be true in the 
general case, as is easily observed if a third constraint is considered together 
with the two others in (4.14) to give: 

f, « -2 ♦ 6 < o 

f 2 = z - 20 - d + 2 < 0 

f 3 = -z ♦ 60 - 9d < 0 (4.16) 

1 < 6 < 2 

The feasible region for this set of constraints is shown in Figure 4.3a 
and Figure 4.3b for d = 1.0, and the corresponding function f is shown in 
Figure 4.3c. Note that y(6,d) is no longer monotonic in 6, and is in fact non- 
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Figure 4.3: Feasible region and y(d,fl) for constraints (4.16) with d=1.0 
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diff erentiable at 8 = 9/5, and exhibits two local maxima at 8 = 1 and 0 = 2. 
It is clear from Figure 4.3a and Figure 4.3b that the size (depth) of the feasible 
region decreases at both extreme points, 8 = 1 and 8 = 2, and gets enlarged 
towards the interior point 8 = 9/5. There are in this case two critical points 
to be considered for design, which are in fact the two extreme points of the 
range 1 ^ 8 < 2. This observation about the location of critical points can be 
generalized for the case of a set of nonlinear convex constraint functions 
through the following theorem: 

Theorem 2. If the constraint functions fj(d,z,0) 

are jointly convex in z and 8, then the problem 



max min max f.(d,z,#) 
0£T z jGJ J 



(4.17) 



has its global solution 8 C at an extreme point of the 
polyhedral region T = { 8 \ 8 L < 8 < 8 U ). 



Proof: 



This theorem can be proved in three parts as follows: 



Property 1. If for every j £ J, fj(d,z,#) is jointly convex in z and 8, 

then <f>(ti,z,8) = max f.(d,z,0) is also jointly convex in z and 8. 
jGJ J 

Proof: Since for any ti, fj(d,z,0) is jointly convex in z and 8, the epigraph 



Epi { f.(d,z f 0)} = { y.z.0 | y£f.(d,z,0)} 



is a convex set, for every j G J 
Hence, 



P| Epi | f.(d,z,0) | is also a convex set 
(see e.g. Stoer and Witzgall, 1970; Rockefellar, 1970). 
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But, P) Epi { fj(d,z,^) } = { y.z.6 | y ^ f Ad.z.6) , V j € J } 
j€J 

= { y,z,0 | y > max f .(d,z,0) } 

= Epi | max f.(d,z,#) j 
j€J 

= Epi { ?(d,z,8) } 

which is therefore a convex set. From this it follows that the function f>(d,z,6) 
is also jointly convex in z and 0. 

Property 2. If </>(c\,z,d) is jointly convex in z and 0, 

then f{d,d) = min ^(d,z,0) is convex in 6. 
z 



Proof: Let jHd,0 3 ) = min $»(d,z,<9 3 ) = j>(d,z 3 ,0 3 ). 



12 3 
Let 6 , 0 £ T be two distinct points that are different from 6 , 

and 0 < X < 1, such that 0 3 = (1-X) 6 1 + X 6> 2 . 

Let ^(d,6> 1 ) = min j6(d,z,6> 1 ) = j»(d,z 1 ,i9 1 ), 
z 

and j/<(d,0 2 ) = min f{6.z,d 2 ) = <f>(d,z 2 ,$ 2 ). 
z 

Since ^(d,z,0) is jointly convex in z and 

(l-x) ^(d,z 1 ,0 1 ) ♦ X j*(d,z 2 ,0 2 ) > ,*(d.z l2 ,0 3 ), 

where z 12 = (1-X) z 1 + Xz 2 . 

But, j>(d,z 12 ,0 3 > > min j>(d,z ( 0 3 ) = j>(d,z 3 ,0 3 ), 



and therefore, (l-x) ^d.z 1 ,0 1 ) ♦ X f(d,z 2 ,6 2 ) > ^(d,z 3 0 3 ), 
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or. jMd,* 1 ) + X f(6,d 2 ) > f(6,d\ 

Noting that (9 3 = (1-X) 0 1 + X 8 2 , it is clear that f(6,6) is convex in 0. 

Property 3. If f(d,6) is convex in 6, 

then every local solution 6° of the problem 

max f(d,d) 

lies at an extreme point of the convex region T, 
unless the solution is degenerate. 

Proof: Assume that 6° is a non-extreme point of the region T. 

1 2 

Then it is possible to choose two distinct points 6 , 0 € T in the 
neighborhood of 0° , and 0 < X < 1 such that 6° = (1-X) 0 1 ♦ X <9 2 . 
Since f(d,d) is convex in 6, 



(l-x) fid.dh + X yid.e 2 ) > fid.e°). 



Since 6° is a non-degenerate solution for the above maximization problem, 



max f(6,6) > (l-x) jMd,0 1 ) ♦ X f(6,6 2 ) 



Therefore, max f(6,d) > f(6,d°), 



which is a contradiction since 6° maximizes locally the function <f>(d,d). 
Hence, the assumption that d° is a non-extreme point of the region T must be 
incorrect, and this proves the result stated above. 
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Property 4. If the region T is polyhedron defined by (4.1), 

the global solution 6 C for the problem (4.17) 

must lie at a corner point (vertex) of this polyhedron, 

unless the solution is degenerate. 

This result is obvious from the fact that the vertices are the only extreme 
points for a polyhedral region, and that every boundary point (as well as 
interior point) can be expressed as a convex combination of these extreme 
points (vertices). Therefore, any local solution to the problem (4.18), and 
hence its global solution d c , must lie at a vertex of the polyhedral region 
T. Thus the result stated in Theorem 2 is proved. 

4.5 DISCUSSION 

Since there are a finite number of vertices for the poiyhedron T, 
Theorem 2 provides an answer to the question as to whether a finite number 
of points can be considered for design to ensure feasibility for all the points 
in the polyhedron T. It follows from Theorem 2 that if the constraints are 

convex, feasibility of operation for every value of 6 6 T can be guaranteed by 
considering in the design all the vertices of the polyhedron T, since all the 
critical parameter values d c will then be accounted for. Also, since (4.17) 
represents a maximization of a convex function as shown in Property 3 of 
Theorem 2, there can be multiplicity of local solutions for (4.17), and hence a 
number of different critical points. This result indicates that in optimizing a 
design, it is important to ensure feasibility for every one of the possibly 
several critical points. The solution algorithms presented in the next chapter 
provide a systematic procedure to ensure feasibility of the design for every 
critical parameter point by explicitly considering them for design. 

It should also be clear that the assumption of convexity on the 
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constraint functions in Theorem 2 is a sufficient condition for the location of 
the critical points at the vertices of the polyhedron T. Therefore, there can 
also be cases where even if non-convex constraint functions are involved, the 
critical points correspond to vertices. However, it is clear that this will not 
always necessarily be true. Despite this limitation, the result provided by 
Theorem 2 is very useful for deriving efficient algorithms for solving the 
problem of design under uncertainty, as shown in the next chapter. 
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CHAPTER 5 



SOLUTION ALGORITHMS 
FOR DESIGN UNDER UNCERTAINTY 



As shown in Chapter 4, when solving the problem of optimal design 
under uncertainty, it is essential to satisfy the max-min-max constraint in order 

to guarantee feasibility of operation for every 6 G T. This implies that if for a 
given design d, feasibility holds for every critical parameter value d c , then 

feasibility can be guaranteed for every SET. In this chapter, two solution 
algorithms are presented which assume that the critical points are located at 
the vertices of the polyhedral region T of the parameters. These algorithms 
transform the problems of design under uncertainty into an iterative multiperiod 
design problem, for which an efficient decomposition strategy has already 
been presented in Chapter 2. Finally, a brief discussion is presented on the 
problems and limitations of locating critical parameter points. 
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5.1 ALGORITHM 1. 



As was proved in Theorem 2, if the inequality constraint functions f. , 
j € J are convex, then the critical points must lie at the vertices of the 
polyhedron T. Since there are a finite number of vertices in T, a design 
obtained by considering all these vertices will be feasible for every point in 
the polyhedron. This would then suggest the following algorithm for solving 
the design problem: 

Step 1 - Include all the N vertices in the set 

T Q = { | e x is a vertex of T, i=1,2 n] (5.1) 

Step 2 - Solve the problem 

N 

minimize C°(d) + jT w l C(d,z',0') 
d. z\z 2 z N ft 

s.t. f(d,z',0') < 0 i = 1,2 N (5.2) 

e T o 



with the projection-restriction strategy, 
so as to obtain the design d°. 



Since T Q includes all . the N vertices of the polyhedron T, every critical point 
d c corresponding to the above design d° will also be included in T . 
Therefore the design d° will be feasible for its critical points, and hence it 
will be feasible for every 6 E T. 

The drawback in this algorithm is that the number of vertices N to be 
considered for design increases exponentially with the number of parameters p, 
since N = 2 P . Thus, for a problem involving ten uncertain parameters (p = 10), 
the design problem has to consider 2^ - 1024 vertices, which would lead to 
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an extremely large problem in (5.2). Because of this fact, Algorithm 1 is only 
suitable when the number of parameters is small (typically, not more than 
five). 

5.2 ALGORITHM 2. 

When the number of uncertain parameters is large, it is desirable to have 
an efficient solution procedure which does not have to consider all the 
vertices explicitly as in Algorithm 1, but that can still ensure the optimality 
and feasibility of the design. The following algorithm is proposed for this 
purpose, wherein only a subset of the vertices are considered for design in 
the resulting multiperiod formulation. 
Step 1 - Set k = 0. 

Choose an initial set T Q consisting of N Q vertices where N Q « 2 P . 

This can be achieved with small computing requirements, using the 
procedure suggested by Grossmann and Sargent (1978), in which each 
constraint is maximized individually by assuming monotonicity and analyzing 

the signs of their gradients. The gradients dfJdd^ of each of the individual 
constraint functions f , , j=1,2,....m, with respect to the parameters 6^ ,k=1,2,....p, 
are computed at initial values of d and z, and the signs of these gradients are 

analyzed. If for each individual constraint function f, , the gradient 

dfJdQfc > 0, the upper bound 0^ ' s selected for the parameter 6^, whereas if 
dfj/d0j< < 0 the lower bound 9^ is selected. Clearly, for zero gradients either 
choice of the bounds is possible. Since each constraint may lead to a 
different vertex, the set of vertices obtained for all constraints is finally 
merged into into the smaller set of vertices T Q by using a set covering 
formulation. It should be noted that // the constraint functions f ; , are 
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monotonic in the parameters 6^, these vertices will correspond to the 
maximization of individual constraint functions (see Grossmann, 1977). 
Step 2 - Solve the problem 



N, 



minimize C°(d) + y w 1 m,z\e') (5.3) 
d, z',i=1,2 N k 1-1 

s.t. f(d,z',0') < 0, i = 1,2 N k 



so as to obtain the design d . 

c k 

Step 3 - Determine the critical parameter values 8 ' by solving for every 
vertex #' in T^, the problem 



y(d k ,0 ] ) = min{ u | u ^ f .(d k ,z,d ] ), j Gj} (5.4) 
z 

The vertex that gives rise to the maximum value of f is then 
determined and is denoted by 0 c ' k . If f(6 k ,6 z,k ) < 0, stop. 
Otherwise, proceed to Step 4. 

Step 4 - Define T R+1 = T k U { 6 cM } , N k+1 = | T k+1 |, (5.5) set 

k = k+1 and iterate from Step 2. 

». • 

Note that at the termination of this algorithm the design will necessarily 
be feasible for all values of parameters, because it will be feasible for the 
critical parameters values. Also, the algorithm has to terminate in a finite 
number of iterations since there are only a finite number of critical parameter 
points to be considered. The initial vertices predicted in Step 1 by the 
method of Grossmann and Sargent (1978) will often yield very good guesses 
for which only one iteration in Algorithm 2 may be required. As in 
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Algorithm 1, problem (5.3) in Step 2 can be solved with the projection- 
restriction strategy described in Chapter 2. Also, note that the minimizations 
in (5.4) may not have to be performed until completion for all vertices, as 
they can be stopped when f reaches a negative value in which case the 
existence of a non-empty feasible region is detected. Thus, by the above 
considerations Algorithm 2 will provide in general a much more efficient 
method of solution than Algorithm 1. 

However, there are two factors in Algorithm 2 which would require 
further investigation. One of them is the number of parameter points 
considered for design, which in turn determines the size of problem (5.3) in 
Step 2. This number will increase at each iteration since a new parameter 
point will be added. The question is whether this number can be kept small 
throughout, by eliminating some of the previous points while adding new ones. 
This elimination can probably be performed on the basis of the value of f, to 
remove those vertices corresponding to the minimum value of f. From the 
discussions in Chapter 4 it is clear that a minimum value of f indicates a 
feasible region with the largest depth and hence most likely to remain feasible 
for small changes in the design d. It should be noted that no theoretical 
justification is given for the elimination of vertices in this manner, and it is 
to be considered only as a very good heuristic. A second question is whether 
it is possible to determine the critical parameter points in Step 3 without 
explicitly analyzing each of the individual vertices and solving for f(ti,d). If 
this were possible, it would further enhance the efficiency of the above 
solution strategy when dealing with a large number of parameters. 
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5.3 NUMERICAL EXAMPLES 

To illustrate the application of the algorithm described above, two 
example problems are presented below. 



In this example the heat exchanger network 4SP1 of Lee et al. (1970) with 
outlet temperatures specified as inequalities is considered (see Grossmann and 
Sargent, 1978). The flowsheet consists of five heat exchangers, one of which 
is a steam heater and another being a cooler using cooling water, and with 
two hot streams and two cold streams as shown in Figure 5.1. Table 5.1 
gives the data for the problem. The overall heat transfer coefficients Uj, 
i=1,2,...5, were considered to be the parameters with ±20% uncertainty in their 
values. The design problem then consists in selecting the areas Aj, i=1,2,....5 
so that irrespective of the actual values of the heat transfer coefficients 
(within the ±20% range around their nominal values), the specifications on the 
outlet temperatures are satisfied by suitable choice of the cooling water outlet 
temperature T^ and the steam temperature T^. Apart from the equality 
constraints representing the heat balance and design equations for the network, 
the following inequality constraints on the temperatures of various streams 
have to be satisfied. 



Example 1. 
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T 11 " T 5 * °' 55 

T y - T 2 > 0.55 
T g - T 6 > 0.55 
T g - T 15 > 0.55 

Among the above constraints, the first four constraints correspond to 
specifications on the outlet temperatures, and the last five on the minimum 
temperature approach. Table 5.2 gives the initial set of vertices considered 
for design, which were obtained by analyzing the signs of gradients of 
individual constraints as suggested by Grossmann and Sargent (1978). Note 
that this set includes also the nominal point so as to provide an adequate 
weighting in the cost function. The CPU-time required to solve this 
multiperiod design problem with the projection-restriction strategy was only 
16.7sec (DEC-20). The design corresponding to these five parameter points 
was found to be feasible for all the 32 vertices. The results are given in 
Table 5.3, from which it is clear that the values of f are non-positive at all 
the vertices, thus ensuring feasibility. 

Note that the actual value of f is dependent on the scaling factors used 
for the constraint functions, which are given in Table 5.1 for this problem. 
Although the choice of these scaling factors is arbitrary, it does not affect the 
detection of (in)f easibi I ity for a given design. The CPU-time to determine f 
for the 32 vertices and thus test for feasibility was only I8.0sec (DEC-20). It 
is also interesting to note that there were many vertices with the same value 
of f(d,d), which may be attributed to the fact that the constraints themselves 
are sparse with respect to these uncertain parameters. In solving this example 
problem, optimization was performed in each case with the variable metric 
projection algorithm of Sargent and Murtagh (1973). In this example the 
optimum design solution was obtained with a single iteration through 



- KpH - 



- 78 - 



Algorithm 2, without the need for considering additional vertex points. This 
may not always be the case as can be seen from the second example which 
follows. 
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Table 5.1: Data for Example 1. 

Nominal values of the uncertain parameters: 

U 1 = U 2 = U 3 = U 5 = 3066 kJ / m 2 hr K 

U 4 = 4088 kJ / m 2 hr K 

Other process parameters: 

FC 1 = 21912 kJ / hr K T 1 = 389 K 

FC 4 = 27461 kJ / hr K T 4 = 333 K 

FC y = 38009 kJ / hr K T ? = 522 K 

FC n = 31674 kJ / hr K T n = 433 K 

T 14 = 311 K 

Cost function: 

5 5 

C = 145.6 ^A°- 6 + r w 1 ( 18.5 F l s * 0.923 F 

i=1 i=1 

o» 1 = 0.6, a- 1 = 0.1. i=2,3,4,5. 
Scaling factors for constraints: 1.8 

Bounds on control variables: 

314 ^ T 15 < 355 K 

534 < T 13 < 556 K 
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Table 5.2: Parameter values considered for design in Example 1. 



u 1 u 2 u 3 u 4 u 5 

1 N N N N N 

2 U U U U N 

3 L U U U N 

4 U L L L L 

5 L U L U N 



N - Nominal, L - Lower bound U - Upper bound 
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Table 5.3: Results for Example 1. 
2 

(a) Heat exchange areas, m 

A 1 A 2 A 3 A 4 A 5 cost, $/yr 

30.8 62.2 45.58 3.9 2.9 1 1758 

CPU-Time (DEC-20) for obtaining the design 
using the projection-restriction strategy: 16.7sec 

(b) Test for feasibility at the vertices: 



f(6,6) Number of vertices Vertex number, v 
0.0 24 0-15, 24-31 

-0.679 Nominal Point 

-2.452 4 20-23 

-3.775 4 16-19 



CPU-Time (DEC-20) for testing feasibility 

and determining f(d,d) at the vertices: 18.0sec 



Vertex number v given by: 

v . T r, (2) 5H . 
i=1 



•{ 



o, if u, = ir 



1. if Uj = U 



u 
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Example 2. 

Figure 5.2 shows the flowsheet consisting of a reactor and a heat 

exchanger, used to handle a first-order exothermic reaction A B, for which 
the problem data is given in Table 5.4. The parameters considered to have 
uncertainty in their values are: (i) F Q , the feed flow rate (±10%); (ii) T , the 
temperature of the feed stream (±2%); (iii) T +, the inlet temperature of 
cooling water (±3%); (iv) k R , the Arrhenius rate constant (±10%); and (v) U, the 
overall heat transfer coefficient for the cooler(± 10%). Among these five 
parameters, the first three are associated with inlet streams to the units while 
the latter two correspond to internal parameters of the process. The 
conversion is specified to be not less than 90%, and the temperature of the 
reactor must be lower than the specified upper bound, 389°K. The design 
problem consists in selecting the optimal sizes for the reactor and heat 
exchanger so that the specifications can be satisfied by suitable choice of the 
temperatures T2, 7^2' in F '9 ure ^.2, irrespective of the actual values of 
the parameters. The material and heat balance equations and design equations 
for the reactor and heat exchanger represent the equality constraints of the 
design problem, and are similar to the equations (3.9H3.14). Other 
specifications to be satisfied are expressed by the following inequality 
constraints: 

V 

0.90 

389 
0 

0 



(a) 
(b) 

(c) 
(d) 



V > 



( C Ao " C Al) ' 



Ao 
333 < 



T 1 " T 2 * 



(e) 



T w2 " T w1 ~ 
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(f) 



T, - T. 



w2 



> 11.1 



<g> 




W1 



> 11.1 



The initial set of parameter points consists of the nominal point and 
three vertices, obtained by analyzing the gradient of the constraints, as given 
in Table 5.5a. The design corresponding to these four points is given in 
Table 5.6a. This design is found to be infeasible for eight of the thirty-two 



vertices. Since the value of f is found to be the same for all these eight 
vertices, one among them is chosen to be added to the initial set of vertex 
points considered for design. This new set of parameter points is given in 
Table 5.5b and the resulting design shown in Table 5.7a. This design is 
feasible for all the 32 vertices as shown by the non-positive values of y 
given in Table 5.7b. Here again, these values of f correspond to the scaling 
factors given in Table 5.4 for the constraints of the problem. This example 
illustrates the need to analyze the max-min-max constraint as a means to 
achieve feasibility of operation for the specified range of parameter values. 
Again, as in the previous example, there are many vertices having the same 
value of f(d,d), which can be attributed to the fact that the constraints are 
sparse with respect to the parameters. Note that it took 9.2sec and 12.8sec to 
solve the multiperiod v design problem with 4 and 5 periods respectively, in 
Step 2 of Algorithm 2. As can be seen from Tables 5.6 and 5.7, the 
computational time taken to test for feasibility was comparatively large - 
65.8sec and 73.9sec respectively in the two iterations through Algorithm 2. It 
took a total of 161.7 sec (DEC-20) for obtaining the complete solution, which 
is quite moderate, considering the size and complexity of this problem. The 
optimizations were performed using the variable metric projection algorithm of 
Sargent and Murtagh (1973). 



vertices as indicated in Table 5.6b by the positive values of f for these eight 
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Finally, the sensitivity of the solution with respect to different choices 
of weighting factors (in the objective function) was tested for this problem. 
Table 5.8 gives the results obtained for this design problem by taking different 
weighting factors w l in the objective function. As can be seen from this 
table, the volume of the reactor is not dependent on the weighting factors, 
whereas the heat transfer area and the annual cost varies within ±11.5% and 
±5.7% respectively. It is to be noted that the design does not get 
considerably affected for moderate changes in the weighting factors. It is 

only when aP -» 1.0 that the variations seem to be quite substantial, since in 
this case most of the weighting is given to the nominal point. Although the 
actual design solution depends on these weighting factors, the feasibility of 
the design is always guaranteed. This is because the constraints 
corresponding to the critical parameter points are satisfied at the solution, and 
it is only through the sensitivity of the objective function that the weighting 
factors affect the optimum design solution. However, this sensitivity is 
usually not large since the design has been optimized for several parameter 
points, unlike the common practice which is to optimize for a single parameter 
value. 
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Figure 5.2: Reactor - Heat Exchanger System for Example 2. 
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Table 5.4: Data for Example 2 



Nominal values of uncertain parameters: 

k R = 0.6242 hr" 1 U = 1635 kJ/m 2 hr K 

F Q = 45.36 kgmole/hr T Q = 333 K 
* T w1 = 300 K 

Other process parameters: 

E/R = 555.6 K " AH rxn = 23260 kJ/kgmole 

C Aq = 32.04 kgmole/m 3 C = 167.4 kJ/kgmole 



Cost function: 



N 

C = 691.2 V 0 - 7 ♦ 873.6 A 0,6 + «' ( 1.76 + 7.056 F\ 

i=1 

« 1 = a; 0 ; « (1 - w°)/(N-l), i=2,3 N 

uP - 0.5, (weighting factor for the nominal point) 

Scaling factors for constraints: 

(a) 3.531; (b). 100.0; (c) 1.8; (d),<e).(f).(g) 18.0 

Bounds on control variables Starting Point 

V = 7.09 m 3 

A = 11.15 m 2 

311 < T 1 < 389 K T 1 367 K 

31 1 < J 2 < 389 K 7 2 = 333 K 

301 - T w2- 355 K T w2 = 350 K 
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Table 5.5: Parameter values considered for design in Example 2. 



k R U F o T o T w1 



N N N N N 

L U L L L 

U L U U U 

U L U U L 

L L U U U 



(a) Initial set of points: (1), (2), (3), (4) 

(b) Second set of points: (1), (2), (3), (4), (5) 



N - Nominal, L - Lower bound, U - Upper bound 
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Table 5.6: Results for Example 2, First iteration with Algorithm 2. 

(a) Design obtained for parameter points given in Table 5.5a 

V = 5.3 m 3 , A = 10.5 m 2 Cost: 10820 $/yr 

CPU-Time (DEC-20) for obtaining the design 
using the projection-restriction strategy: 9.2sec 

(b) Test for feasibility at the vertices: 



f{6,9) Number of vertices Vertex number, v 

+ 1.280 8 4-7, 12-15 

0.0 Nominal Point 

0.0 16 0-3, 8-1 1, 20-23, 28-31 

-1.151 8 16-19, 24-27 



CPU-Time (DEC-20) for testing feasibility 

and determining f{d,6) at the vertices: 65.8sec 



Vertex number v given by: 

v = y. j (2) 5 -» , <r={ 

1=T i. if e. = 

^ = k R , * 2 = u, * 3 = F 0 , * 4 = T o< * 5 et w1 
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Table 5.7: Results for Example 2, Second iteration with Algorithm 2. 



(a) Design obtained for parameter points given in Table 5 




3 9*66. /eftoz 

V = 6.5 rrT, A=«ft*-m z Cost: - ¥) I W $/yr 

CPU-Time (DEC-20) for obtaining the design using 
the projection-restriction strategy: 12.8sec 

(b) Test for feasibility at the vertices: 



f(6,6) Number of vertices Vertex number, v 

0.0 8 4-7, 12-15 
-1.220 Nominal Point 

-1.220 16 0-3, 8-11, 20-23, 28-31 

-2.323 8 16-19, 24-27 



CPU-Time (DEC-20) for testing feasibility and 
determining f(d,8) at the vertices: 73.9sec 



Vertex number v given by: 
v = TV (2) 5 -' , 

0, = k R , 8 2 = U, ^ 3 = F 0 , 0 4 = T Q , T w1 
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Table 5.8: Results of Example 2, for different choices of weighting factors 



0 

w 


V m 3 


A m 2 


Annual Cost, $/yr 


0.0 


6.5 


10.22 


1 1396 










0.1 


6.5 


10.13 


11302 


0.2 


6.5 


10.02 


1 1206 


0.3 


6.5 


9.91 


1 1 108 


0.4 


6.5 


9.79 


1 1007 


0.5 


6.5 


9.66 


10902 


0.6 


6.5 


9.51 


10793 


0.7 


6.5 


9.35 


10678 


0.8 


6.5 


9.14 


10553 


0.9 


6.5 


8.86 


10409 


1.0 


6.5 


8.11 


10160 



Weighting factors for the vertices: 

1 0 
w = w ; (Nominal Point) 

J = (1 - a/°)/(N-D, i=2,3 N 

Number of periods considered N = 5. 
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5.4 DISCUSSION ON LOCATING THE CRITICAL PARAMETER POINTS 

It is interesting to see that in the first example the design solution is 
obtained by considering only those vertices that were indicated by analyzing 
the gradients of individual constraint functions. As pointed out earlier in 
Section 5.2, if the constraint functions are monotonic with respect to the 
parameters, this procedure corresponds to the maximization of the individual 
constraint functions (see Grossmann and Sargent, 1978). However, it is clearly 
evident from Example 2 above that such procedure cannot always lead to the 
right set of parameters points required to be considered for design. The only 
method of characterizing these parameter points is by analyzing the max-min- 
max constraint. Thus, Example 2 illustrates this need for analyzing the max- 
min-max constraint as a means to achieve feasibility of operation for the 
specified range of parameter values. To understand this aspect in the overall 
design strategy, consider the following set of three linear constraints: 

f 1 = -z 1 + 3(9 1 - e 2 < 0 

f 2 = -z 2 - B y ♦ 30 2 < 0 

f 3 = z 1 + z 2 - 6 y - 6 2 - d < 0 (5.6) 
1 < 0j < 2, i=1,2. 

Since these constraint functions are linear, their gradients are independent 
of the initial point chosen for such calculations; and since they are monotonic 
in d, the maximization of each of these functions can be performed by 
analyzing these gradients. It is clear that the three constraint functions get 
maximized at the three different vertices 0 1 = [2,1], 0 2 =[1,2] and 0 3 =[1,1] 
respectively. Therefore, according to the procedure which is based solely on 
maximization of individual constraint functions, these three vertices would 
form the set of parameter points to be considered for design. The fact that 
this procedure can lead to designs that may be infeasible is shown below by 
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e 2 = [1.2] 




Figure 5.3: Feasible Region for the constraints in (5.6) for d = 3. 
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analyzing the max-min-max constraint for the above set of design constraints, 
which determines the value of f - the measure of (in)feasibility for any given 
design d. For a design d, the value of f(d,9) for any parameter value 6 is 
given by 

f(d,d) - min u 
z 

s.t. f- = -z- ♦ 30- - 6 2 < u 

f 2 = ~ z 2 " d -\ + 26 2 ~ u (5 ' 7) 

f~ = z- + z 0 - 0. - #~ - d < u 

Figure 5.3 gives a plot of the feasible region for the set of constraints in 

(5.6), with a design d = 3, wherein the values of f obtained from (5.7) are also 

given for each of the four vertex points of the parameter-space. Since the 

constraint functions are linear and hence convex, the critical parameter points 

must lie at a vertex, as shown in Theorem 4.2. The values of f in Figure 5.3 

therefore indicate that #^=[2,2] is the critical point, since f attains its 

maximum at this vertex. Also note in Figure 5.3 that the design d = 3 is 

1 2 

found to be feasible for the three vertices 6 =[2,1] , 6 =[1,2] and 
6 =[1,1], whereas it is infeasible for the critical parameter point - namely 
the vertex # 4 =[2,2]. To make the design feasible for the critical point as 
well, consider d = 4. The feasible region and the value of f for the four 
vertices are again shown in Figure 5.4 for this value of d = 4. It is clear from 
Figure 5.4 that the design is just feasible for the critical point, whereas a 
finite region of feasibility exists for all other points in the parameter space. 
Thus, by ensuring the feasibility of the design for the critical parameter 
values, it is possible to guarantee the overall feasibility of the design for 
every parameter value within the specified range. 

In order to gain some further insight as to why in some cases the 
maximization of individual constraint functions will lead to the correct critical 



- KpH - 



- 96 - 



points, assume that for all 0ET the same common (single) set of values z 

is selected for the control variables. It then follows that the max-min-max 
constraint reduces to 

max max f.(d,z,#) ^ 0 

del jEJ J 

which is equivalent to (5.8) 

max fAd.ieh < 0, y j E J. 

Thus if for the design design d, it is possible to select a control z, feasible 

and common for all 9 E T, then some of the parameter points predicted by 
maximization of individual constraints will correspond to those (critical points) 
predicted by the max-min-max constraint. However, it is clear that in general 
different controls z may have to be selected for different realizations of 8 to 
maintain feasibility. Therefore, maximization of the individual constraint 
functions cannot always lead to the critical points, although it usually provides 
a very good guess for an initial set of parameter points to be considered for 
design. 

Another interesting question about locating the critical points is to 
determine the conditions under which only a single critical point would suffice 
to be considered for design. As pointed out earlier in Chapter 4, monotonicity 
of f(6,d) in 0 gives rise to a single critical point in the space of parameters, 
as is also the case in the example given above in (5.6). Although the property 
of monotonicity of f is difficult to establish when the constraint functions 
are nonlinear, the following analysis can be applied for the case of linear 
constraint functions. 
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Suppose that for a given design d, the set of linear constraint functions 
are expressed as: 

P n z 

v2I a jk*k +2>jk z k +c j * °- i- 1 - 2 m (5 - 9) 

k«1 k=1 
The critical point 0 is given by the solution to the problem: 



max min u (5.10) 

eei z 

P n z 

s.t. f . - X a jk^k + X b jk z k + c j " u < i" 1 ' 2 m 

k=1 k=1 



Assume that for any 6 £ T the problem 

f(6,6) = minu (5.11) 

2 

P n z 

sa - f j = X a jk*k + X b jk z k + c j - u ' j= 1 - 2 m 

k=1 k=1 

has the same set of active constraints (e.g. the first r, r ^ m). 

The Kuhn-Tucker conditions for the above minimization problem then 

gives 



(a) 1 = T 

j=1 



r 

X. 



r n z 



(b) 0 = X. b jk , k = 1,2 n 2 (5.12) 

j=1 k=1 



(c) 



P 

I . 

k=1 k=1 



- Z a jk *k + H b jk £ k + c j ■ i* 1 - 2 m 
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Since the value of u at the minimum determines f, it follows that 

r P n z 

f - TH(Z a ik*k + Z b jk z \ +c j) < 5 - i3) 

i=1 k=1 k=l 

From this expression it is clear that f(d,d) is monotonic in 6, for the chosen 
design d. Note that in (5.12) the active constraints are determined by the 
values of the multipliers \y which are in turn obtained from the subset (a),(b) 
in (5.12). Since there are r multipliers X., and n 2 + 1 equations in (5.12a, b), 
and since r ^ m, a necessary condition for having the same set of 

constraints to be active for every 8 E T would be: m ^ n z + 1. In general, a 
necessary and sufficient condition would be to have the values of all the 
multipliers X. (j = 1,2,....m) uniquely determined by the system (5.12). 

Furthermore, if for every k, s\gn(dy / ldd ^) ■ sign(a^) for some constraints j, 
then these constraints are maximized at the point defined by max f{d,d), and 
under these conditions the maximization of individual constraints with respect 
to 6 would lead to the critical points. 

To illustrate these ideas consider the set of constraints given by (5.6). 
The solution to (5.12 a,b) yields X^ = \^ = X^ = +1/3, indicating that the 
three constraints are active for all d. Furthermore, from (5.13) the value for f 
is obtained as f{ti,d) = 1/3 ( d - + 6 2 " d ) » indicating that fid, 8) is indeed 
monotonic in 8, and that maximization of f results in the critical point given 

by # 4 = [2,2]. Again, the signs of the gradients df/d8^ and djW8#2 are both 
positive, whereas there is no constraint in (5.6) that satisfies this condition. 
Therefore, in general there is no way of relating the monotonicity of f with 
the monotonicity of the individual constraint functions f. . From this it 
should also be clear why the maximization of individual constraint functions 
does not predict the right critical parameter values. It is also important to 
note here, that in general the set of active constraints would depend on the 
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design d, from which it follows that the monotonicity of f(6,6) may very well 
depend on the design d, when m > n 2 + 1. For instance, in the set of 
constraints (4.16), f is not monotonic in 6 for d = 1, as shown in 
Figure 4.3c. However, as can be seen from Figure 5.5, for d = 0.5, f is 
monotonic in 6. 

Unfortunately, the analysis presented above would seem to be too 
restrictive to be applicable for the general case of nonlinear constraint 
functions. However, the analysis does provide some insight on the condition 
for monotonicity of f{d f 9), which might suggest that it would be worth to 
investigate further this point. 
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5.5 NOMENCLATURE FOR THE NUMERICAL EXAMPLE PROBLEMS 



Heat transfer area of the heat exchanger, m 
Annual cost, $/yr 

3 

Concentration of reactant in the feed stream, kgmole/m 

Concentration of reactant in the product stream, kgmole/m' 
Heat capacity of reaction mixture, kJ/kgmole K 



Ratio of activation energy to gas constant, K 

Feed flowrate, kgmole/hr 

Flowrate of recycle, kgmole/hr 

Flowrate of steam, kg/hr 

Flowrate of cooling water, kg/hr 

Heat capacity flowrate, kJ/hr K 

Heat of reaction, kJ/kgmole 

Reaction rate constant, hr" 1 

Total number of vertices considered for design, (integer) 
Temperature of feed stream, K 
Reaction temperature, K 
Recycle temperature, K 
Temperature of stream i, K 
Temperature of cooling water, K 
(subscript 1 - inlet, subscript 2 - outlet) 



Overall heat transfer coefficient, kJ/m hr K 

3 

Reactor volume (design capacity), m 
3 

Reaction volume, m 

w ] : weighting factor corresponding to vertex i. 

Superscript i refers to the vertex i (period of operation), i=1,2,....N 
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CHAPTER 6 

CONCLUSIONS 
AND RECOMMENDATION FOR FURTHER RESEARCH 



The problem of optimal design of flexible chemical plants has been 
studied in this thesis. In order to account for changes in operating conditions 
and yet meet the design specifications, systematic procedures have been 
proposed to introduce the required flexibility in the optimal design of a 
chemical plant. It is shown that there can be two classes of problems in 
optimal design of flexible chemical plants - the deterministic multiperiod 
problem and the problem of design under uncertainty. The problem 
formulation and solution strategies for both these classes of problems have 
been discussed. 



- KpH - 



- 103 - 



6.1 DETERMINISTIC MULTIPERIOD PROBLEM 

If the operating conditions are given as a discrete sequence, the problem 
of optimal design can be formulated by using the deterministic multi period 
model. This leads to a nonlinear program with block-diagonal structure in the 
constraints, wherein the main computational difficulty arises because of a 
large number of decision variables involved. It is shown that existing 
decomposition techniques cannot be applied effectively to circumvent this 
problem. A new decomposition scheme has been proposed here, based on a 
projection-restriction strategy which exploits very effectively the mathematical 
structure of the problem and the fact that many inequality constraints become 
active at an optimum design solution. The performance of this projection- 
restriction strategy has been analyzed to show that significant gains in 
computational effort can be achieved with this technique. The relationship 
between the performance of the projection-restriction strategy and the size 
(number of periods) of the design problem being solved is discussed, wherein 
the effect of having many active constraints at the solution is clearly shown. 
The results of the numerical example show a linear increase in computational 
time with the number of periods. It is interesting to note that the proposed 
decomposition scheme does not presuppose the use of any particular type of 
optimization algorithm for solving the nonlinear programming problem, and is 
conceptually quite general in its application for the type of design problems 
discussed. 

As is shown, the efficacy of the proposed decomposition technique is 
dependent on the number of inequality constraints that actually become active 
at the optimum design solution. It should be noted that this is not a serious 
limitation of the technique, since it is a common observation in design that 
many inequality constraints do in fact become active at the optimum solution 
of a design problem. On the other hand, the decomposition technique exploits 
this additional feature of the design problem, thus resulting in a very efficient 
solution strategy. 
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It is mentioned that there are two important questions to be considered 
for a successful implementation of the proposed projection-restriction strategy. 
One of them is the problem of finding an initial feasible point for the 
problem. The proposed method for this purpose consists in replacing the 
objective function in the original design problem by the sum of squares of the 
constraint violations, and solving it by optimizing alternatively with respect to 
the design and control variables. This is found to be very efficient since it 
does decompose the problem to a manageable size, and handles an objective 
function that is relatively well behaved. 

The second question is concerned with a procedure for eliminating the 
variables in the restriction step, by using the active constraints, for reducing 
the number of decision variables of the nonlinear program. This involves the 
selection of an appropriate set of new control variables so that the resulting 
system of equations is solvable. Here again, a procedure has been proposed, 
which is in fact shown to be applicable also in the general context of selection 
of decision variables and detection of redundancy in rectangular systems. 

6.2 DESIGN UNDER UNCERTAINTY 

When the values of some of the parameters are either known only 
approximately or they actually vary or fluctuate during the operation, the 
design problem has to be formulated so as to account for this uncertainty in 
such parameters. Since it is usually difficult to obtain a prior knowledge of 
the probability distribution for these parameters, before the implementation of 
a design, it is assumed that the range of values for these parameters is 
available. However, in cases where the probability distribution is known, the 
bounds on these parameter values may be obtained as appropriate confidence 
limits for these values which can be determined from such data. 

Having fixed the range of values for these parameters, which are 
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assumed to be independent, the objective is to propose a design strategy that 
can guarantee feasible steady state operations of the plant for every such 
value of these parameters. That is, the plant is designed with the required 
flexibility, so that depending on the actual values of these parameters (within 
the given bounds} the operation of the plant can be suitably adjusted to meet 
the design specifications. A rigorous mathematical formulation which 
represents this design strategy has been shown to be a two-stage nonlinear 
infinite program, wherein optimization is to be performed on the set of design 
and control variables so as to satisfy the constraints for every possible value 
of the parameters within the specified bounded polyhedral region, while 
minimizing the expected value of the annual cost function. To guarantee 
feasibility of the design for every possible realization of the parameters, a 
logical constraint is incorporated in the problem. This feasibility constraint 
ensures that for every allowable value of the parameters that may be 
encountered during the operation, appropriate values of control variables can 
be chosen so as to satisfy the design constraints. 

The main computational difficulty with the above formulation arises 
because of the infinite number of variables and an infinite number of 
constraints involved. By selecting a discrete set of parameter values and 
approximating the expected value of the cost function by a weighted average, 
the number of control variables can be reduced to be finite. To circumvent 
the problem of an infinite number of constraints, an equivalence for the 
feasibility constraint has been established, which leads to a max-min-max 
constraint. This max-min-max constraint has also proved to be very useful in 
obtaining a deeper understanding of the feasibility aspect of the design 
problem. The following are the main results obtained from an analysis of the 
properties of the max-min-max constraint in relation to the design problem. 

• It is shown that for a given design and a fixed parameter value, the 
max-min-max problem provides a measure of (in)feasibility of the 
operation. 
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• A critical parameter value corresponds to the smallest degree of 
feasibility for operation, or the largest degree of inf easibility. 

• In general there can be more than a single critical parameter value 
that must be considered for design, in order to ensure feasibility. 

• If the constraints are convex the critical points must lie at the 
vertices of the polyhedral region of parameters. 

Based' on these results an algorithm has been proposed which transforms the 
problem of design under uncertainty into an iterative multiperiod design 
problem, wherein each period would correspond to operating the plant with a 
distinct value of the parameters. The proposed algorithm leads to an efficient 
solution procedure, as has been shown by two example problems. 

6.3 RECOMMENDATIONS FOR FURTHER RESEARCH 

In the formulation of the deterministic multi period design problem the 
constraints have been considered to be coupled in the design variables, but 
uncoupled in the state and control variables among the different periods. If in 
fact there is a subset of constraints interconnecting these variables in 
different periods, the problem will then have a bordered block-diagonal 
structure in the constraints. However, the proposed projection-restriction 
strategy could still be applicable if in the linking constraints only a few 
coupled state and control variables occur, as then these can be treated as 
design variables to recover the block-diagonal structure. This problem of 
coupling constraints needs further work so as to devise a scheme which can 
render the projection-restriction strategy applicable, while still maintaining its 
efficiency as a promising solution procedure. Additional work is required also 
for a theoretical investigation of the convergence properties of the proposed 
projection-restriction strategy. It should be noted here, that Ritter (1973) has 
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presented a proof for his version of the decomposition technique when applied 
to linear constraints. Regarding further numerical experience with the 
projection-restriction strategy, a general computer program incorporating this 
strategy has been the objective of a Master's thesis project (Avidan, 1982). 
Such a capability will certainly be useful for further investigations into the 
possible refinements of the proposed strategy so as to make it applicable 
also for the more general case of multiperiod design problems with coupling 
constraints. 

In the problem of optimal design with uncertain parameter values, one of 
the factors in the proposed solution Algorithm 2 that requires further 
investigation is the question of avoiding the cumulative build-up of the number 
of parameters to be considered for the multiperiod design problem in Step 2 
of Algorithm 2. It is suggested that this could be achieved by removing at 
every iteration of Algorithm 2, at least one parameter point, on the basis of 
the minimum value of f{d,d), which indicates a parameter value with a feasible 
region having the largest depth and hence most likely to remain feasible for 
small changes in design. This procedure would seem to be a very good 
heuristic, in the absence of a more efficient and theoretically valid approach 
which can still guarantee feasibility for the entire range of parameter values. 
Further investigation is needed in this direction, to study the scope and 
limitations of this approach. 

As for Step 3 of Algorithm 2, it would be desirable to devise a more 
efficient procedure for solving the max-min-max problem when a large 
number of parameters is involved. Note that in the proposed algorithm the 
value of f{6,6) has to be determined explicitly for every vertex point in the 
given polyhedral region of the parameter space, in order to find the vertex 
that maximizes f{6,9). Therefore, this approach for locating the critical points 
could become quite inefficient as the number of uncertain parameters 
increases, since the corresponding number of vertices increases exponentially. 
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One possible solution strategy would be to linearize the constraints and 
perform a parametric analysis of the resulting linear program to determine the 
critical parameter points, along with an appropriate procedure for bounding the 
value of f to account for the error in the above linearization. However, to 
make this a valid procedure for finding critical points in the nonlinear case 
would seem to be a nontrivial task which would require considerable research 
work. 

Note that the uncertain parameters are considered to be varying 
independently of one another, both in the mathematical formulation of the 
design problem in Chapter 4, as well as in presenting the solution algorithms 
in Chapter 5. It should be pointed out here, that the two-stage programming 
formulation given in Chapter 4 is in fact quite general, and is valid even when 
the parameter variations are linearly constrained (through their interdependence). 
Also, the basic idea of the proposed solution algorithm can indeed be 
applicable for the case of interdependent parameter variations, if a suitable 
procedure can be developed for locating the critical points in Step 3 of 
Algorithm 2. However, this would involve the ability to identify the extreme 
points of the constrained convex region of parameters, and in general it would 
be desirable to locate the critical points through a complete solution for the 
max-min-max problem. 

It is also important to note that the assumption of convexity of the 
constraint functions is only a sufficient condition for the location of critical 
points to be at the vertices (extreme points) of the given polyhedral (convex) 
region in the parameter space. In fact, a weaker condition would be to relax 
the assumption of convexity of the constraint functions, and instead impose 
the condition of convexity of the feasible region. Also, since in the 
nonconvex case the critical points may not correspond to vertices, it would be 
desirable to devise a procedure for solving explicitly the max-min-max 
constraint. However, the complexity of this constraint, and the fact that it can 
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lead to a nondifferentiable problem with multiple local solutions would appear 
to make it very difficult to have an efficient solution procedure valid for the 
general case. 
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6.4 CONCLUSIONS 



This research work has resulted in the following contributions to 
the area of Computer-Aided Chemical Process Design: 

*• 

A systematic approach for introducing flexibility in optimal design 
of chemical plants, that is more rational than other existing 
approaches. 



• A very efficient decomposition scheme based on a projection- 
restriction strategy for solving multiperiod design problems, which 
results in significant gains in computational effort. 

• A rigorous mathematical formulation for the problem of optimal 
design with uncertainty in parameter values, for which its 
interpretation and properties have been investigated. 

• An efficient strategy for solving the problem of design under 
uncertainty that results in designs which can be guaranteed to be 
both feasible and optimal. 

• A procedure for detection of redundancy and selection of decision 
variables in process design computations, so as to avoid the 
problems of singularity, in solving systems of equations using 
tearing procedure. 
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APPENDIX 

SELECTION OF DECISION AND TORN VARIABLES 
IN PROCESS DESIGN COMPUTATIONS 



One of the main requirements of chemical process computations is the 
solution of sparse systems of nonlinear algebraic equations. For instance, it 
has been mentioned in Chapter 2 that the state variables are determined by 
solving the set of equations representing the chemical process system, which 
are in general nonlinear. Also, in Step 3 of the projection-restriction 
strategy, the active constraints are added to the set of equality constraints so 
as to define a new system of equations, for which a selection of decision 
variables (control variables) has to be obtained. In this case, it is in fact 
possible that there are more equations than there are variables, and therefore 
it is desirable that the procedure being developed can also handle such 
situations. However, very often these systems have several degrees of 
freedom since the number of variables is usually larger than the number of 
/ independent equations. In order to derive a solution procedure for these 
rectangular systems, it is necessary to identify redundant equations in order to 
determine the number of degrees of freedom of the system, and to select a 
suitable set of decision variables, so that the resulting square system of 
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equations is solvable. This Appendix presents such a procedure which is 
applicable not only to the projection-restriction strategy, but also to general 
design computations. 

A.1 INTRODUCTION 

In formulating a solution procedure for systems of algebraic equations 
that arise in process design computations, the set of variables is partitioned 
into a vector of decision variables and a vector of state variables. The state 
variables are calculated by solving the system of equations for any fixed 
values of decision variables. These values are chosen either by the engineer 
or through an optimization of some appropriate objective function. 

It should be noted that the determination of the number of degrees of 
freedom requires the identification of a set of consistent and non-redundant 
equations. This in general is a non-trivial problem since the number of 
equations for a chemical process system can be rather large. Therefore, it is 
important that automated solution procedures perform this task so as to allow 
the engineer to specify all the valid and relevant equations, even if some of 
these equations are redundant. 

A common approach for deriving a solution procedure for these 
rectangular systems of equations consists in ordering the equations, selecting 
decision variables and tearing state variables, so that the occurrence matrix of 
the system has a bordered lower triangular form, as shown in Figure A.I. In 
solving the system of equations for any fixed values of decision variables, the 
torn variables are iterated while the remaining state variables are calculated by 
solving the equations in the lower-triangular portion. Clearly, it is desirable to 
select an ordering of equations so as to minimize the number of torn 
variables. A number of algorithms based on this approach have been reported 
in the literature. See, for example, Lee et al., (1966); Lee and Rudd, (1966); 
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Figure A.1: Bordered lower-triangular structure of the system (A.2) 
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Christensen and fiudri, (1969); Christensen, (1970); Westerberg and Edie, (1971); 
Leigh, (1973); Stadtherr et al., (1974); Hernandez and Sargent, (1979). When the 
system has redundant equations, or if there are singularities associated with 
certain choices of decision variables, all these algorithms can fail to obtain a 
square system which is non-singular. In fact, only Edie and Westerberg (1971) 
present an algorithm for choosing the decision variables to avoid singularity of 
the resulting square system, but they assume that the original system does not 
contain redundant equations. The method described below is suitable for 
detecting redundancy in a given set of equations and selecting an appropriate 
set of decision variables, so that the resulting system is solvable, without the 
problem of inconsistency and singularity. 

A.2 CONSISTENCY AND SINGULARITY 

For the purpose of the present discussion, let the system of nonlinear 
algebraic equations be 

h(x) = 0, x G X C R n , f c R m (A.1) 

In order to analyze the system of equations it is useful to define the three 
properties: structural (in)consistency, algebraic (non-)singularity and numerical 
(non-)singularity. The system is said to be structurally consistent if by 
permutations of the rows and columns of its occurrence matrix it is possible 
to arrive at a rearrangement of the occurrence matrix having non-zero entries 
along the diagonal. In other words, the system is structurally consistent if 
and only if it has at ieast one admissible output set, as shown by Steward 
(1962). Structural inconsistency arises because of the presence of a subsystem 
having more equations than the number of state variables defining the 
subsystem, and therefore always results in a jacobian matrix that is rank 
deficient. 
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While I inlconsistency is a structural property of a system, (non-)singularity 
is associated with the algebraic and numerical characteristics of a system. 
The system is algebraical ly singular if and only if its jacobian matrix is rank 
deficient at all points in the domain X of the system. If the jacobian matrix 
is rank deficient only at certain points being analyzed, then the system is said 
to be numerical ly singular at those points. Note that algebraic singularity 
does not necessarily imply structural inconsistency. An example of this is a 
consistent system of equations which are linearly dependent. Also note that 
algebraic non-singularity is indicated by the presence of at least one point 
where the system is numerically non-singular. Thus, for any system of 
algebraic equations the following implications hold: 

structural ly inconsistent =^ algebraical ly singular 

=^> numerically singular at every point in the domain X. 

and, 

numerically non-singular at every point in the domain X 
z£> algebraically non-singular structurally consistent. 

From these implications, it is clear that structural consistency is only 
necessary but not sufficient to ensure algebraic non-singularity. For successful 
numerical solution, most algorithms (e.g. Quasi-Newton methods) require that 
the system of equations be at least algebraically non-singular in the state 
variables in the region of computation. Algebraic singularity may be caused 
by structural inconsistency or by the presence of redundant equations or even 
by a bad choice of decision variables. Thus structural analysis alone is not 
sufficient for deriving a procedure to eliminate redundant equations and to 
select a right set of decision variables. As shown below, it is possible to 
modify existing equation ordering algorithms so as to ensure algebraic non- 
singularity in the system of equations. 
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A.3 IDENTIFYING REDUNDANCY AND SELECTING DECISION VARIABLES 

First, the rows and columns of the occurrence matrix of the equations in 
(A.1) are permuted using any of the above cited algorithms, so as to obtain a 
bordered lower triangular form. The reordered system of equations f'(x') = 0 
can be represented as shown in Figure A.2, by the two sets of equations: 

s(u,v) = 0 

r(u,v) = 0 (A.2) 

where [f'] T = [s T , r T ] and [x'] T = [u T ,v T ]. 

Here, s is the set of non-recycle equations with lower triangular structure 
in the state variables v, which can therefore be solved sequentially for the 
vector v given a value of u, and r corresponds to the set of recycle 
equations. Since v can be treated as an implicit function of u, the system of 
equations in (A.2) can be reduced to the form: 

r(u,v(u)) = 0 (A.3) 

The vector u includes the set of decision variables d as well as the vector of 
torn variables t, that is, u"^" = (d" 1 "^" 1 "). For any chosen value of d, the system 
may be solved by iterating on the torn variables t, to reduce the residuals of 
the equations in (A.3) .to zero. The partitioning of the vector u is in general 
not trivial, and requires careful consideration to avoid singularities, as shown 
below. 

The subsystem s in (A.2) can be assumed to be algebraically non-singular 
because of its lower-triangular structure. Hence, the algebraic singularity of 
the system can be analyzed through the jacobian J c (r,u) of (A.3), which is in 
fact the constrained jacobian of the original system (A.2) corresponding to the 
recycle equations. This matrix is given by: 
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Figure A.2: Constrained jacobian matrix for the system (A.2) 
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J c< r ' u > ' ¥ - ¥ <A.4) 

c 3u dv N dv ' du 

and is indicated in Figure A. 2 by the entire shaded area. In order to analyze 
the singularity of the system in (A. 3), we need to determine the rank of this 
jacobian matrix J (r,u). With this, a corresponding non-singular submatrix J* of 
highest rank can be found as shown in Figure A.2 by the darker area. The 
variables in u that form the columns of this submatrix J* will then be chosen 
as the torn variables t. The remaining variables in u will correspond to the set 
of decision variables d, and its dimensionality determines the number of 
degrees of freedom of the system. Those equations r r that are not included 

in the rows of this submatrix J* will be deleted, since they are either 
redundant or inconsistent. Thus it can be ensured that the jacobian matrix of 
the resulting sauare system is of full rank and hence non-singular (Halemane 
and Grossmann, 1981). Since there can be several possible choices for the 
submatrix J*, it is also possible to incorporate any preferred choice for the 
decision variables in the above procedure. 

In practice, the jacobian matrix J c <r,u) given in (A.4) can be evaluated 
numerically at the current point by performing small perturbations in the vector 
u. Since equation ordering algorithms tend to minimize the dimension of the 
vector u, the work involved in this procedure is often very moderate. The 
rank of the jacobian matrix J (r,u) and a corresponding non-singular square 
submatrix J* of highest rank can be obtained by performing a Gaussian 
elimination on this matrix. For this elimination, the procedure presented by 
Westerberg et al. (1979) can be applied. Also, in order to ensure that the 
redundant equations are always deleted and the remaining system is non- 
singular, the specified tolerances for the pivots should be larger than the 
maximum round-off error that may be accumulated in the elimination process. 
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It is to be noted here that by analyzing the jacobian (A.4) at a single 
point, there is sometimes the danger of not being able to distinguish between 
a non-singular jacobian matrix that is ill-conditioned and a truly singular matrix. 
This can cause the deletion of a non-redundant equation and introduce more 
design variables than required. To guard against this possibility, one may 
evaluate the jacobian at different points. However, a better strategy would be 
the following procedure: 

Whenever a numerical solution is obtained for the square system, the 
residuals of the deleted equations are evaluated to check whether they are 
satisfied within a given tolerance. If any of these deleted equations is not 
satisfied, it indicates that the number of degrees of freedom of the system 
has been over-estimated, and therefore requires some decision variables to be 
converted to state variables. This can be done by a reanalysis of the 
constrained jacobian at the current point. It should be noted that one can 
expect in general to perform only one analysis of the constrained jacobian, 
except for severely ill-conditioned problems. 

A.4 EXAMPLE 

To illustrate the ideas presented above, the example of a binary 
isothermal flash system (Edie and Westerberg, 1971) is considered, as shown 
in Figure A.3. This can be represented by the set of equations (a) to (I) as: 



P° - P°(T) 
P° - P°(T) 



= 0 



= 0 



(b) 



(a) 




= 0 



= 0 



(d) 



(c> 




= 0 



(e) 
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^2 


- 


^2^2 




■ 0 


(f) 


Z 1 


* 


z 2 " 


1 


: 0 


(g) 


x 1 


+ 


x 2 - 


1 


' 0 


(h) 


v 1 


- 


y 2 " 


1 « 


■ 0 


(i) 




+ 


y,v - 




1 0 


(j) 




+ 


y 2 v - 


z 2 F . 


' 0 


(k) 


L 


+ 


V 


F 


0 


(i) 



In this case all the equations are clearly valid, but they contain one redundant 
equation, for instance, equation (I). However, this equation will be included in 
the system so as to illustrate the proposed procedure presented earlier. 

Table A.1 shows a structure obtained by equation ordering, with the 
entries corresponding to the non-zero elements of the jacobian matrix of the 
original system. The constrained jacobian J Q (r,u) corresponding to the recycle 
equations, as given by the expression in (A. 4) is obtained algebraically and 
shown in Table A.2. It is clear that the rank of this matrix is two, as it can 
be seen in Table A.2 that the row (k) is obtained by adding row (I) to L times 
row (h). Hence, one of these three recycle equations is redundant and needs 
to be deleted. Table A. 3 gives the possible choices for decision variables, 
that arise by the deJetion of any one of these three equations. Further 
analysis as for a final selection of a set of decision variables can be made 
depending upon the requirements of the particular problem at hand. 
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A.5 DISCUSSION 

The above example shows how the proposed method can be used to 
identify redundant equations and to select a proper set of decision variables, 
so as to result in a non-singular square system which can be solved 
numerically using suitable equation solving algorithms. Also, it is important to 
note that, by this method one is able to determine the right number of 
degrees of freedom of a system, which is essential in both the analysis and 
design of any process system. In using the proposed method, the engineer is 
allowed to specify all relevant and valid equations that describe the system, 
even if some of the equations are redundant. The analysis required to identify 
redundancy and to select a right set of decision variables as presented in the 
proposed method, can be incorporated in an automated solution procedure. 

The numerical information on the jacobian matrix, which is required for 
the analysis presented earlier, can in fact be used also as an initial estimate 
in the updating formula of a Quasi-Newton method (e.g. see Broyden, 1965) 
when solving the recycle equations. The proposed method in itself, is not 
limited by the numerical evaluation and analysis of the jacobian matrix at any 
given point. This is clear from the example wherein no numerical evaluation 
has been used. If an automatic algebraic manipulating package is assumed to 
be available, it should be possible to detect the algebraic singularity of the 
jacobian matrix, and hence of the system. It should be noted here that 
although this procedure seems to be directly applicable only when using a 
tearing procedure for equation solving, there is no reason why it cannot be 
extended for other types of equation solving algorithms. 
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Figure A.3: Binary isothermal flash system 
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Table A.1: Structure of the binary isothermal flash system 



as obtained by equation ordering 



P° P° K 1 K 2 y 2 x 1 x 2 z 2 L V F z 1 y 1 P T 

(a) 1 a 

(b) 1 b 

(c) 1 -P -K 1 

(d) 1 -P -K 2 

(i) 1 1 

(e) -x 1 -K 1 1 

(f) -x 2 1 -K 2 

(g) 1 i 

(j) L x 1 Yl -z 1 -F V 

(h) 1 1 

(k) V L -F x 2 y 2 -z 2 

(I) 1 1 -1 
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Table A.2: Constrained Jacobian J c <r / u) 
for the structure given in Table A.1 

_>v 

V F 2l P T 



(h) 0 0 0 1/K 1 - 1/K 2 1/P ax^K^ + bx 2 /K 2 P 

(k) 1 - y 1 /x 1 -1+z 1 /x 1 F/x 1 -(V + L/KgJ/x^ 0 -Lax^H^P + Lbx 2 /K 2 P 
(I) l-y 1 /x 1 -1+Z 1 /X 1 F/x 1 -(V + L/K 1 )/x 1 -UP - La/K^P 
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Table A.3: Possible choices of decision variables for 
the ordering given in Table A.1. 



Equation 


Decision 


Torn 


deleted 


Variables 


Variables 


(h) 


V. F. z v y, 


P, T 




V, F, z y P 


y r t 




V F 7 P 
» » r, l y r 






V, F, P, Yl 


z r T 




V, F, T, 


Zy P 




V, F, P, T 






V, Zy Yy P 


F, T 




V, Zy Yy T 


F, P 




V, Zy P, T 


f, y t 




F, Zy Yy P 


V, T 




F. Zy Yy T 


V, P 




F, z r P, T 


V, y, 


(k) or (I) 


V, F, Zy y 1 


P. T 




V.*F. z r P 


y r t 




V, F. Zy T 


Yy P 
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22 MONDAY 

TARTAN. The Tartan meets every Monday 
evening at 9 p.m. in Skibo room 13. Anyone 
interested in writing, photography, or working 
on any other facet of newspaper production 
is welcome to stop by or call the offices at 
x21 12. 

SUPA SEMINAR Dr. Keith Poole will speak on 
"Congessional Voting" at 3:30 p.m. in MMC 
room B-14. 

JOINT PHYSICS COLLOQUIUM. Professor Frank 
von Hippie of Princeton University will speak 
on "Physicists in the Public Policy Debate", at 
4:30 p.m. in WEH 7500. Coffee will be 
served at 4:15 p.m. Everyone is welcome. 

23 TUESDAY 

JOINT LOGIC SEMINAR. CMU and the 
University of Pittsburgh will sponsor a seminar 
at 2 p.m. in WEH 5336. For more 
information, check the main bulletin board 
outside WEH 6114. 

SOCIAL SCIENCE BROWNBAG. Steve Klepper, 
Steve Garber and Dan Rubenson will discuss 
"The Underground Economy" at 12:30 p.m. in 
PH 217. Everyone is welcome. 

METALLURGICAL ENGINEERING AND 

MATERIALS SCIENCE GRADUATE SEMINAR. 
The seminar will feature J. Gilbert and will be 
held at 3:30 p.m. in DH A3 17. Abstracts will 
be located on the bulletin board outside WEH 
3325. 

STUDENT SENATE. The Student Senate meets 
at 7 p.m. in the Wherrett Room, Skibo. 
Everyone is welcome and students may 
address the senate during audience 
participation. 

NUMERICAL ANALYSIS WORKSHOP. A 
numerical analysis workshop will be held in 
WEH 6121 at 9 am. For more information, 
see the bulletin board outside WEH 6113. 




MATHEMATICS GRADUATE SEMINAR. A 
graduate seminar will be held at 5:30 p.m. in 
WEH 5310. See the bulletin board outside 
WEH 61 13- for details. 

A PHI O. Alpha Phi Omega, the national service 
fraternity, meets every Tuesday evening at 
6:30 p.m. in Skibo 12. All are welcome. 

24 WEDNESDAY 

MECHANICAL ENGINEERING SEMINAR. J.M. 
Beer of MIT will discuss "Fluidized Combustion 
of Coal" at 2:30 p.m. in Scaife Hall Auditorium. 
Everyone is welcome to attend. 

FICTION READING. Bruce Dobler will present a 
public reading of his short fiction at 8 p.m. in 
Baker Hall 231. Everyone is invited to attend. 

PERSPECTIVES ON WORLD CITIES SERIES. 
"Urban Economics" will be discussed at a 
lecture by Professor William Alonzo of 
Harvard University from 4-5 p.m. in Gregg 



Hall. The lecture is free and open to the 
public 

ART LECTURE SERIES. George Trakas, a 
sculptor from New York City, will present his 
personal views on art at 8 p.m. in Gregg Hall. 

PH.D THESIS DEFENSE. K.P. Halemane will 
defend his thesis, "Studies in the Optimal 
Design of Flexible Chemical Plants", at 12:30 
p.m. in the Chemical Engineering Conference 
Room Open to the Public. 

EPP SEMINAR. V Penninti will speak on 
"Fundamental Studies on Hot Cathode Arcs" in 
Scaife Hall 301 at 12 noon (Bring a lunch.) 

AIESEC-CMU. The International Association of 
Students in Econqmics and Business 
Management welcomes students of every 
major to their weekly meeting at 4:30 p.m. in 
Skibo 1 1. 

THREE STOOGES FILM FESTIVAL Kappa Sigma 
fraternity will sponsor a Three Stooges Film 
Festival in DH 2210 at 6,8.10,12 p.m. 
Admission is $1.00. 

25 THURSDAY 

MENTAL HEALTH POLICY COLLOQUIUM SERIES 
Dr. Herbert J. Schelsinger, professor of 
psychiatry, School of Medicine, University of 
Colorado, will present a lecture entitled 
"Mental Health Services and Medical Care 
Utilization in the Fee for Service System" at 4 
p.m. in the Multi- Purpose Room of Hunt 
Library. The lecture is free and open to the 
public. Coffee will be served at 3:30 p.m. 

HUMAN SEXUALITY SEMINAR. Dr. Tom Allen, 
gynecologist will discuss "Sex - How to Keep 
it Safe" at 8 p.m. in the Wherrett Room. 
Refreshments will be served. For more 
information, contact Anne Witchner X2107. 

ASTRONOMY CLUB. The Astronomy Club 
meets every Thursday at 7:30 p.m. in DH 
2304. All are welcome. 

26 FRIDAY 

EPP SEMINAR. Erick Erickson will speak on 
"Predator and Prey Relationships in Ecology 
and Economics" at 3 p.m. in SH 206. 

27 SATURDAY 

NO SCHEDULED EVENTS. 

28 SUNDAY 

NO SCHEDULED EVENTS. 



1 MONDAY 

NO SCHEDULED EVENTS. 

2 TUESDAY 

METALLURGICAL ENGINEERING AND 

MATERIALS SCIENCE SEMINAR The seminar 
will feature G. Pierini and will be presented in 
DH A3 17 at 3:30 p.m. Everyone is. welcome 
to attend. 

STUDENT SENATE. The Student Senate meets 
at 7 p.m. in the Wherrett Room, Skibo. 
Everyone is welcome and students may 
address the senate during audience 
participation. 

NUMERICAL ANALYSIS WORKSHOP. A 
numerical analysis workshop will be held in 
WEH 6121 at 9 a.m. For more information, 
see the bulletin board outside WEH 6113. 

MATHEMATICS GRADUATE SEMINAR. A 
graduate seminar will be held at 5:30 p.m. in 
WEH 5310. See the bulletin board outside 
WEH 6113 for details. 

JOINT LOGIC SEMINAR. CMU and the 
University of Pittsburgh will sponsor a seminar 
at 2 p.m in WEH 5336. For more 
information, check the main bulletin board 
outside of WEH 6114. 

A PHI O. Alpha Phi Omega, the national service 
fraternity, meets every Tuesday evening at 
6:30 p.m. in Skibo 12. All are welcome. 

3 WEDNESDAY 

MECHANICAL ENGINEERING SEMINAR. 
D. Anderson of Purdue University will speak 
on "Research in CAD/CAM and Computer 
Graphics" at 3:30 p.m. in Scaife Hall 
Auditorium. Everyone is invited to attend. 

EPP SEMINAR. Mark Shelton will speak on 
"Optimizing Toxicity, Energy Conservation and 
Profitablity in Industrial Chemical Complexes" 
at 12 noon in SH 301 (bring a lunch.) 

ART LECTURE SERIES. Nancy Hagin, a painter 
from New York City, will discuss her personal 
views on art at 8 p.m. in Gregg Hall. The 
seminar is open to the public. 

HIGH ENERGY PHYSICS SEMINAR. Dr. Goran 
Senjanovic of Brookhaven National 
Laboratory will discuss "An Oasis in the Desert: 
Weakly Broken Parity in Grand Unification 
Theories" at 4:15 p.m. in WEH 7500. Coffee 
will be served at 4 p.m. All are welcome. 

AIESEC-CMU. The International Association of 
Students in Economics and Business 
Management welcomes students of every 



major to their weekly meeting at 4:30 p.m. in 
Skibo 11. 

4 THURSDAY 

METALLURGICAL ENGINEERING AND 

MATERIALS SCIENCE STAFF SEMINAR. Dr. 
Gordon H. Geiger of the Research Labs, Intend 
Steel Co. will present a seminar entitled "The 
Thermal Stability of Countercurrent Packed- 
Bed Reactors" at 3:30 p.m. in DH A3 17. For 
more information contact Marge Cain or HI 
Aaronson at x2709. 

HUMAN SEXUALITY SEMINAR. Debbie Smith, 
Sex Consultant will discuss "What Happens 
During Sex?" at 8 p.m. in the Highlander Room. 
Refreshments will be served. For more 
information contact Anne Witchner at x2107. 

JOINT NUCLEAR PHYSICS SEMINAR. Dr. Donald 
Senjanovic of the Argonne National Laboratory 
will speak on "Macroplastic Features of Pion 
Inelastic Scattering and Charge Reactions" at 
4:30 p.m. in DH 1212. Coffee will be served 
at 4:15 p.m. in WEH 8327 (8th floor lounge). 

ASTRONOMY CLUB. The astronomy club meets 
every Thursday in DH 2304 at 7:30 p.m. All 
are welcome. 

5 FRIDAY 

PROFESSIONAL ETHICS LECTURE SERIES 
Douglas Maclean, research associate at the 
University of Maryland's Center for Philosophy 

6 Public Policy, will present a lecture entitled 
"Social Values and Benefit-Cost Analysis" at 
3:30 p.m. in the Wherrett Room. 

SPIRIT FASHION SHOW. "Controversy, a 
Fashion Extravaganza" will be presented in the 
Skibo Ballroom. Dinner will be served at 6:30 
p.m. in the Faculty Dining Room followed at 
8:15 p.m. by the fashion show. Tickets for 
the dinner and show are $7.50, for the show 
only, $4.00. For ticket information call 682- 
9209, 241-6920, or x3924. 

H&SS SAC PARTY. The "Tropics Party" will 
dispel those midsemester chills from 4:30- 
6:30 on the Baker Hall Rotunda. Spring Fling 
tickets will be on sale. All H&SS 
students, faculty, and staff are welcome. 

6 SATURDAY 

NO SCHEDULED EVENTS 

7 SUNDAY 

NO SCHEDULED EVENTS. 



MOVIES 

GERMAN FILM SERIES 

Feb. 23 Una Braake 

7 p.m. 

Feb. 2 Der starke Ferdinand 
7 p.m. 

All films will be shown in Gregg Hall. 
Admission is free, and the interested public is 
welcome. 

THURSDAY FILM ARTS. 

Feb. 25 Breaker M or ant 

6, 8, 10 p.m. 
Mar. 4 Swept Away 

6, 8:30, 1 1 p.m. 

All films will be shown in DH 2210. 
Admission is $1 with ID; $1.25 without. 

AB FRIDAY FILMS 

Feb. 26 Excalibur 

6, 8:30, 1 1 p.m. 
Mar. 5 Arthur 

6,8,10,12 

All films will be shown in DH 2210, 
Admission is $1 with ID; $1.25 without. 

SDC FILMS. 

Feb. 27 A Night at the Opera 

6, 8, 10 p.m. and 12 midnight. 
Mar. 6 S.O.B. 

6, 8:30, 1 1 p.m. 
All films will be shown in DH 2210 at 6, 8, 
10 p.m. and 12 midnight Admission is $1 
with ID; $1.25 without 

SUNDAY EVENING AT THE MOVIES. 

Feb. 28 North by Northwest 
The film will be shown at 6, 8:30, and 1 1 
p.m. in DH 2210. Admission is $1 with ID; 
$1.25 without 

SPORTS 

MEN'S BASKETBALL. 

Feb. 23 Washington & Jefferson 

(A) 8 p.m. 
Feb. 25 Bethany 

(H) 8 p.m. 

WOMEN'S BASKETBALL. 
Feb. 24 Bethany 

(A) 7 p.m. 
Feb. 26 Washington & Jefferson 
(A) 7 p.m. 

SWIMMING. 

Feb. 25-27 

PAC Championships 
Allegheny 



FENCING. 

Feb. 24 California State 

(H) 6 p.m. 
Feb. 27 North Atlantic 

Tournament 

(H) 9 am. 
Mar. 6 Women's Invitational 

IUP 

(A) 9 am. 

GALLERIES AND EXHIBITS 

HUNT INSTITUTE FOR BOTANICAL 
DOCUMENTATION. An exhibit of selected 
works from the recently acquired long-lost 
drawings from the Sesse and Mocino 
expedition to New Spainf 1787- 1803) will be 
on display through March 6. The gallery is 
located in the penthouse, of the Hunt Library 
and is open Monday-Friday, 9 a.m-5 p.m. 

HEWLETT GALLERY. Veteran realist painter 
Malcolm Parcell will exhibit his works from 
Jan. 30 to Feb. 26. The Hewlett gallery is 
located in the College of Fine Arts and is 
open Monday-Friday, 11 a.m.- 5 p.m.; Sat. 10 
a.m -5 p.m.; and Sunday 1-4 p.m. 

HEWLETT GALLERY. Nancy Hagin will exhibit a 
collection of her realistic paintings from 
March 2-21. The Hewlett Gallery is located 
in the College of Fine Arts and is open 
Tuesday to Friday 1 1 am. to 5 p.m., Saturday 
1 1 a.m. to 4 p.m., and Sunday 1 p.m. to 4 
p.m. 

VELAR GALLERY. Sarah Ratchye will exhibit 
her work Feb. 15-26. The Gallery is located 
in the lobby of GSIA and is open Monday- 
Friday, 9 a.m.-5 p.m. 

HUNT LIBRARY FINE AND RARE BOOK ROOMS. 
An exhibition entitled "Recent Gifts" will Vun 
through March 12. The exhibit will be open 
Monday-Friday, 1:30-4:30 p.m. and is located 
on the fourth floor of the Hunt Library 

THEATRE 

MERTON OF THE MOVIES. The world premier 
of "Merton of the Movies" previews Feb. 20- 
24 at 8 p.m., opens Feb. 25-27, March 1-6 
at 8 p.m., with matinee performances on Feb. 
27 and Mar. 6 at 2 p.m. in Kresge Theatre 
The musical version of "Merton of the 
Movies" has a book by Mel Shapiro and Simon 
Burke. The music was composed by Gary 
Friedman, composer of the Broadway 
productions "The Me Nobody Knows" and 
"Platinum". Lyrics are by Bob Lorick. For 
additional ticket information, please call the 
Kresge Theatre Box Office at x2407. 

ON BROADWAY PLAYERS. Scotch N' Soda 
presents "You're A Good Man, Charlie Brown" 



on Mar. 4-7 at 8 p.m. in GSIA Auditorium. 
Tickets are $2.50. 

MUSIC 

FACULTY RECITALS 

Mar. 9 Louis Paul , clarinet 

Ralph Zitterbart, piano 
Arthur Kubey, bassoon 
ACH 8 p.m. 

GRADUATE RECITALS 

Feb. 28 Patricia F/orig, piano 

ACH 3 p.m. 
Mar. 6 Paula Tuttfe, Cello 

ACH 8 p.m. 

STUDENT RECITALS 

Mar. 7 Craig Johnson, Flute 

ACH 8 p.m. 
Mar. 7 Donald Zimmerman, organ 

Kent Smith, French Horn 

Calvary Episcopal Church 

4 p.m. 

Mar. 10 Matthew Litter ini 

French Horn 

ACH 8 p.m. 
Mar. 1 1 Sharon Sciore, Trumpet 

ACH 8 p.m 

COMING EVENTS 

DISTINGUISHED LECTURE SERIES. Ansley 
J. Coale, demographer, 4 p.m. Mar. 8, Skibo 
Ballroom. 

CHEMICAL ENGINEERING SEMINAR. Mr. Tom 
Tulig 1 1 am, Mar. 9 DH 1 102. 

H&SS SPRING FLING Friday April 2, 7 p.m., 
Carnegie Music Hall. 

CAMPUS CALENDAR DEAD- 
LINES 

The Campus Calendar will be issued every 
other week. The deadline for entries is ten 
days prior to the Monday of that calendar's 
release and distribution. Submit copy to the 
Campus Calendar, Public Relation's Office, 16u 
Baker Hall, CMU. Below is a listing of the 
deadlines and the Calendar dates 
corresponding to the deadlines for this 
semester. 

Deadline Calendar 
Feb. 26* Mar. 8-28 
Mar 12 Mar 29- Apr. 2 

Apr. 2 Apr. 12-25 

Apr. 16 Apr. 26-May 9 

Apr. 30* May 10-Aug. 29 
(summer calendar) 
*Please note that this calendar covers an 
extended period of time due to recesses 
and/or exams. 
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1 Louis James Hater, Pittsburgh 
B.S.E.E., Carnegie-Mellon 
University, 1976 
M.S.E.E., Carnegie-Mellon 
University, 1978 
Major: Electrical Engineering 
Subject of Dissertation: Data- 
Memory Allocation in the 
Di stributed Logic Desi gn Style 

Keshava Prasad Halemane, 
Mangalore, India 
B.E., Karnataka Regional 
Engineering College, 1972 
M.Tech., Indian Institute of 
Technology, 1973 
M.S., Carnegie-Mellon University, 
1979 

Major: Chemical Engineering 
Subject of Dissertations Studies in i 
the Optimal Design of Flexible 
irrn i n i i p i u ii Im 

Marc Robert Halpern, Franklin 
Square, N.Y. 

Sc.B., Brown University, 1975 
M.S., Carnegie-Mellon University, 
1977 

Major: Civil Engineering 
Subject of Dissertation: Dynamic 
Response of a Saturated 
Poroelastic Halfspace To Imposed 
Surface Loadings 

Neal Michael Holtz, Ottawa, Ontario, 
Canada 

B.Sc, University of Alberta, 1969 

M.Eng., Nova Scotia Technical 

College, 1970 

Major: Civil Engineering 

Subject of Dissertation: Symbolic 

Manipulation of Design 

Constraints 

Ronald Lasser, Buffalo, N.Y. 

B.S., Carnegie-Mellon University, 
1975 

M.S., Carnegie-Mellon University, 
1977 

Major: Mechanical Engineering 
Subject of Dissertation: The 
Effects of Secondary Flow in a 
Converging Channel of Arbitrary 
Curvature and in the Wake 
Downstream of a Cascade of 
Airfoils ^^^^^^^^ 0 _ mm ^^ 

*Kincho Henry Law, Hong Kong 
B.S., B.A., University of Hawaii, 
1976 

M.S., Carnegie-Mellon University, J 
7979 

Major: Civil Engineering 
Subject of Dissertation: Graph- i 
Theoretic Approaches to 
Structural Analysis 



Ying-Kuo Lee, Taipei, Taiwan 
M.S., Massachusetts Institute of 
Technology, 1975 
Major: Mechanical Engineering 
Subject of Dissertation: Shear 
Bands in Elastic-Perfectly Plastic 
Media 

Jumlong Limtragool, Hua-Hin, 
Thailand 

M.S., Carnegie-Mellon University, 
1977 

M.E., Chulalongkorn University, 
1976 

Major: Mechanical Engineering 
Subject of Dissertation: 
Investigation of the Three- 
Dimensional Behavior of a Crack 
Tip Singularity by Holographic 
Interferometry 

•Michael Harvey Locke, Rando' 
Mass. 

B.S., University of Massachusetts, 
1976 yy 
M.S., Carnegie-Mellon University, 
1978 

Major: Chemical Engineering 
Subject of Dissertation: A CAD 
Tool Which Accomodates an 
Evolutionary Strategy in 
Fnginperir^a Dftsjqn Calculations 

* *Luiz Carlos Barreto Martins, 

Bictoriada Conquista Bahia, Brazil 
Met.E., Escola Politecnica da 
U.S. P., 1973 

M.S., Carnegie-Mellon University, 
1975 

Major: Metallurgical Engineering 
and Materials Science 
Subject of Dissertation: Kinetics of 
Reduction of Iron Ores by Mixed 
Gases Under Pressure 

* Stephen Roy McConnel, Emmett, 

Idaho 

S.B., Massachusetts Institute of 
Technology, 1976 
M.S., Carnegie-Mellon University, 
1978 

Major: Electrical Engineering 
Subject of Dissertation: The 
Implementation and Performance 
of a Fault Tolerant Multiprocessor 

Jack Davenport McDonald, Jr., 
Glenolden 

B.S., Carnegie-Mellon University, 
1973 

M.S., Carnegie-Mellon University, 
1974 

Major: Electrical Engineering 
Subject of Dissertation: Target 
Selective Adaptive Control Model 
for Human Eye Movements 



* Michael Carlton McFarland, Boston. 

Mass. 

A. B., Cornell University, 1969 
M.S., Carnegie-Mellon University, 
1979 

Major: Electrical Engineering 
Subject of Dissertation: 
Mathematical Models for Formal 
Verification in a Design 
Automation System 

* *Harry Gilbert Miller III, Columbia, 

Md. 

B. S.E.E., University of Maryland, 
1974 

M.S., Johns Hopkins University, 
1977 

Major: Engineering and Public 
Policy 

Subject of Dissertation: The 
Reliability and Cost Impacts of 
Alternative Levels of 
Centralization in Electric Power 
Systems 

* *Shantilal Manekchand Mohnot, 

Bombay, India 

B. Tech., Indian Institute of 

Technology, 1973 

M.S., Kansas State University, 

1977 

Major: Chemical Engineering 
Subject of Dissertation: Fast 
Radial Mass Transfer in 
Microporous Flow-Through 
Electrochemical Reaction 

John Charles Molburg, Pittsburgh 
B.S.Eng., University of Illinois, 
1971 

Major: Engineering and Public 
Policy 

Subject of Dissertation: 
Mathematical Models for 
Integrated Environmental/ 
Economic Analysis of the Coal-to- 
Electricity System 

* *Kwoon-Chuen Mui, Kowloon, Hong 

Kong 

B.S., Illinois Institute of 

Technology, 1975 

M.S., Carnegie-Mellon University, 

1978 

Major: Chemical Engineering 
Subject of Dissertation: Stability of 
Displacement Fronts in Porous 
Media — Large Finger Analyses 

* Surajit Mukherjee, Calcutta, India 

B. Tech., Indian Institute of 
Technology 

M.S., Iowa State University, 1977 
Major: Chemical Engineering 



